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Subject:  New  information  on  the  characteristics  of  HF  atmospherics 

radiated  by  lightning  flashes,  when  received  at  over-the- 
horizon  distances  on  a highly  directional  antenna  of  the 
type  used  for  HFDF  and  for  OTH  radar  systems. 


Because  you  routinely  receive  our  reports  on  over-the-horizon  (OTH) 
radar  prepared  under  Contract  N00014-75-C-0930,  we  believe  that  the  en- 
closed technical  report  prepared  under  an  SRI  International  independent 
research  and  development  (IR&D)  project  will  be  of  interest  to  you. 


The  thrust  of  this  work  is  to  understand  better  the  characteristics 
of  HF  electromagnetic  noise  radiated  by  thunderstorms,  as  received  at 
over-the-horizon  distances  by  a highly  directional  antenna.  The  findings 
should  be  of  direct  interest  in  HFDF  and  OTH  radar  applications  whenever 
background  noise  levels  determine  system  performance. 

In  particular,  these  findings  can  affect  the  design  of  nonlinear 
circuits  used  to  reduce  the  effects  of  noise  spikes,  and  the  design  of 
adaptive  antennas  used  to  limit  the  effects  of  spatially  fixed  noise  and 
other  interference.  Here  nonlinear  circuits  include  clippers,  limiters, 
and  noise-excision  techniques  made  practical  by  digital  computers. 

The  report,  "Observed  Characteristics  of  Ionospherically  Propagated 
HF  Atmospherics  from  Normal  and  Severe  Thunderstorms,"  was  prepared  by 
Dr.  Walter  B.  Zavoli,  who  used  the  Wide  Aperture  Research  Facility  (WARF) 
operated  by  SRI  in  support  of  OTH  radar  studies  sponsored  by  the  Office 
of  Naval  Research  and  by  the  Electronic  Systems  Division,  Air  Force  Systems 
Command . 


Some  important  conclusions  of  this  work  can  be  summarized  as  follows: 

• A highly  directional  antenna  can  be  thought  of  as  a different 
way  to  view  the  thunderstorm  environment;  it  has  the  effect  of 
bringing  distant  flashes  close  up.  Waveshape  is  preserved  in 
this  process — at  least  over  bandwidths  of  practical  interest — 
so  the  noise  output  of  a receiver  connected  to  a directional 
antenna  will  be  far  more  impulsive  in  nature,  and  less  Gaussian, 
than  the  output  of  a receiver  connected  to  a less-directional 
antenna. 

• A new  class  of  atmospheric  has  been  identified.  This  class 
seems  to  be  associated  with  the  transition  of  a normal  thunder- 
storm into  a severe  one  (i.e.,  one  which  produces  hail,  funnel 
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clouds,  and/or  tornadoes).  At  close  (line-of-sight)  ranges 
these  atmospherics  tend  to  be  masked  by  other  aspects  of  the 
electromagnetic  radiation  from  storms,  and  thus  seem  to  have 
been  missed  in  prior  studies.  But  at  OTH  distances,  assuming 
sufficient  antenna  directivity,  the  atmospherics  are  distinc- 
tive. Thus  they  offer  the  promise  of  a new  means  for  identify- 
ing, studying,  and  possibly  tracking,  storms  whose  intensity 
can  reach  a point  which  makes  them  a threat  to  lives,  property, 
or  the  success  of  military  operations. 

• A further  unanticipated  finding  is  that  single  lightning  flashes 
from  normal  thunderstorms  cannot  be  considered  to  be  an  electri- 
cal point  source  for  highly  directive  HF  antennas.  For  example, 
radiation  (actually,  a succession  of  discharges)  from  a typical 
lightning  flash  occurring  over  1000  km  from  the  WARF  antenna 
was  shown  to  change  its  bearing  by  roughly  one  degree  during  the 
duration  of  the  flash.  These  bearing  fluctuations  are  presumably 
the  result  of  the  successive  discharging  of  various  charged  vol- 
umes within  a cloud  or  system  of  clouds;  such  behavior  could  have 
a significant  effect  on  the  performance  of  adaptive  antennas  and 
on  the  design  of  algorithms  chosen  to  reduce  or  eliminate  the 
effects  of  spatially  fixed  noise. 

Dr.  Zavoli  will  of  course  welcome  any  comments  or  discussion  of  his 
work.  He  may  be  contacted  by  mail  at  SRI  International,  Building  320A, 
333  Ravenswood  Avenue,  Menlo  Park,  California  94025,  or  by  telephone  at 
(415)  326-6200,  extension  2393. 

Sincerely  yours, 


L.  E.  Sweeney,  Jr. 

Director 
Remote  Measurements  Laboratory 
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ABSTRACT 


\ 

Using  a directional  HF  antenna  with  a nominal  beamwidth  of  0.5^,’ 

A 

atmospherics  from  midwestern  thunderstorms  classified  as  severe,  and 
thus  likely  to  produce  tornadoes,  were  studied  in  California.  The 
antenna  was  scanned  across  regions  within  a one-hop  propagation  distance 
likely  to  develop  tornadoes  at  the  time  of  the  test.  These  regions  were 
selected  on  the  basis  of  National  Weather  Service  severe  weather  fore- 
casts. The  existence  of  strong  propagation  from  these  regions  at  the 
radio  frequencies  used  was  verified  by  means  of  backscatter  soundings. 

Analysis  of  the  data  recorded  in  those  directions  and  periods  for 
which  severe  storm  events  were  reported  led  to  the  discovery  of  a class 
of  radio  atmospherics  significantly  different  from  those  radiated  during 
normal  thunderstorm  activity.  Atmospherics  of  this  new  type  take  the 
form  of  randomly  occurring  impulses  having  certain  distinctive  character- 
istics. A technique  was  devised  to  extract  this  class  of  atmospherics 
from  the  remaining  radio  noise  background.  Based  on  23  hours  of  data, 
taken  on  four  separate  days,  a statistically  significant  correlation 
between  the  reception  of  these  impulses  and  the  presence  of  remote  severe 
thunderstorms  was  found,  x 

Comparison  with  HF  line-of-sight  measurements  performed  by  others 
shows  that,  as  a class,  these  atmospherics  were  very  unlikely  to  have 
been  previously  recognized  as  a significant  subclass  within  the  total 
atmospheric  activity  present  during  severe  storms.  Subsequent  analysis 
found  that  the  temporal  distribution  of  these  atmospherics  can  be  closely 
approximated  by  a Poisson  model  having  a slowly  varying  mean.  It  is  the 
increase  in  this  mean  rate  which  correlates  with  observed  severe  storm 
events  such  as  tornadoes,  funnel  clouds,  and  severe  hail.  Finally,  a 
hypothesis  is  presented  suggesting  concentrated  wind  shears  within  the 
storm  as  the  most  probable  mechanism  for  generating  the  observed  atmo- 
spherics. 


i 


This  research  was  initially  motivated  by  the  observation  of  strong 
and  detailed  atmospheric  waveforms  received  on  the  directional  antenna, 
despite  the  absence  of  local  thunderstorm  activity.  Comparison  of  at- 
mospheric noise  simultaneously  received  with  this  antenna  and  with  an 
omnidirectional  antenna  suggested  that  these  atmospherics  were  from 
normal  thunderstorms  of  remote  origin,  lying  within  the  directional 
antenna's  main  beam.  Direct  verification  was  obtained  when  measurements 
taken  close  to  a New  Mexico  thunderstorm  correlated  with  simultaneous 
measurements  taken  1300  km  away  with  the  directional  antenna. 

The  directional  antenna  was  then  reconfigured  so  as  to  receive  in 
eight  contiguous  directions  simultaneously.  This  revealed  readily 
measureable  changes  in  the  angle  of  arrival  of  atmospheric  noise  radiated 
from  the  many  small  discharges  associated  with  a single  lightning  flash. 
It  would  therefore  appear  that  variations  in  horizontal  location  of  in- 
cloud charge  centers  involved  in  individual  remote  lightning  flashes  can 
be  measured  in  this  way.  The  average  horizontal  extent  of  individual 
flashes  measured  from  one  storm  was  9 ± 4 km. 
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I INTRODUCTION 


A.  Purpose 

The  purpose  of  this  research  was  to  investigate  experimentally 

characteristics  of  high  frequency,  HF,  atmospherics  of  remote  origin. 

A directional  antenna  was  used  to  study  atmospherics  from  severe 

thunderstorms  occurring  several  thousand  kilometers  from  the  receiving 

apparatus.  The  term  "atmospheric"  is  used  here  as  defined  by  Horner — 

that  is,  the  electromagnetic  energy  radiated  from  any  electrical  break- 

1 * 

down  occurring  during  a single  lightning  flash.  A single  atmospheric 
can  be  a simple  impulse,  an  impulse  train,  or  a continuous  waveform. 
Atmospheric  noise  is  the  superposition  of  atmospherics  from  many 
lightning  flashes. 

B.  Motivation 

Individual  atmospherics,  strong  compared  to  the  background  noise 
and  apparently  originating  at  skywave  distances,^  have  often  been  re- 
ceived on  a directional  antenna  situated  in  northern  California.  This 

3 4 

antenna,  part  of  the  Wide  Aperture  Research  Facility  (WARF)  ’ was  de- 
signed by  Sweeney5  and  successfully  used  by  Barnum,6  Washburn7  and 
others  to  investigate  properties  of  the  ionosphere  and  HF  radio  propa- 
gation. For  such  research,  strong  atmospherics  were  a source  of  annoy- 

g 

ance,  and  techniques  such  as  those  described  by  McKinney  and  Zavoli 
have  been  devised  to  mitigate  their  effects.  Considering  the  absence 
of  close-by  thunderstorm  activity,  the  strength  and  detailed  structure 
of  these  atmospherics  (which  were  not  observed  on  an  omnidirectional 
antenna)  suggested  to  the  author  that  the  directional  properties  of  the 
WARF  antenna  might  be  used  to  investigate  remote  atmospheric  electricity. 

* 

References  are  listed  at  the  end  of  this  dissertation. 

^For  the  HF  band,  skywave  propagation  (reflection  from  the  ionosphere) 
can  occur  for  distances  substantially  beyond  the  horizon.  Propagation 
paths  using  a single  ionospheric  reflection  (i.e.,  a one-hop  path)  can 
occur  over  distances  of  200  to  4000  km.  Fundamentals  of  ionospheric 
radio  propagation  can  be  found  in  Davies. ^ 
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In  the  past,  atmospherics  have  been  used  to  study  the  physical 
mechanisms  of  the  lightning  discharge.  These  studies  have  concentrated 
on  line-of-sight  lightning  to  maximize  sensitivity,  eliminate  propaga- 
tion effects,  and  correlate  electrical  measurements  with  optical  records. 
For  the  most  part,  the  component  of  atmospherics  below  100  kHz  has  been 
used  because  of  its  strength  and  the  relatively  small  number  of  impulses 
radiated  per  flash. 

Prior  to  the  current  work,  little  research  had  been  conducted 
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using  remote  (i.e.,  skywave-propagated)  HF  atmospherics.  Horner 
pointed  out  the  difficulties  associated  with  obtaining  meaningful 
records  of  individual  HF  atmospherics  of  remote  origin.  In  general, 
such  signals  are  weaker  than  their  VLF  counterparts.  Since  HF  atmo- 
spherics contain  many  more  impulses  (several  hundred),  the  chances  are 
greater  for  overlap  (the  simultaneous  reception  of  atmospherics  from 
different  flashes).  However,  as  demonstrated  by  the  results  described 
in  this  dissertation  an  antenna  such  as  that  at  WARF  has  sufficient 
directivity  and  sidelobe  rejection  to  overcome  these  difficulties  by 
increasing  signal- to-noise  ratio  (SNR)  and  reducing  event  overlap  on 
in-beam  atmospherics.  These  capabilities  enable  the  WARF  antenna  to 
gather  useful  information  from  remote  lightning  and  thus  study  normal 
and  severe  thunderstorm  activity.  Furthermore,  the  flexibility  to 
point  the  antenna  in  any  direction  of  current  interest  and  take  measure- 
ments from  a distance  avoids  the  difficult  problem  of  trying  to  be  at 
the  right  place  at  the  right  time  (i.e.,  under  an  active  thunderstorm 
of  interest).  This  capability  is  particularly  useful  in  the  study  of 
the  more  elusive  severe  thunderstorms  where  research  is  slowed  by  these 
logistics  difficulties.  The  thrust  of  the  research  described  here  was 
first  a general  investigation  of  the  properties  of  atmospheric  noise 
received  on  the  directional  WARF  antenna,  and  then  a study  of  remote 
normal  and  severe  thunderstorm  activity  through  reception  of  atmo- 
spherics . 
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c. 


Background 


By  the  start  of  this  century,  the  state  of  knowledge  of  atmo- 
spherics paralleled  knowledge  of  the  physics  of  lightning;  each  con- 
tributing to  the  understanding  of  the  other.  Physics,  meteorology, 
powerline  engineering,  and  radio  science  were  the  principal  disciplines 
involved  in  the  study  of  lightning  and  its  electromagnetic  radiation. 

A brief  history  of  early  lightning  research  is  presented  below. 

The  study  of  atmospherics  is  then  surveyed,  with  emphasis  placed  on  HF 
research  pertinent  to  the  current  work.  Finally,  the  previous  research 
on  atmospherics  associated  with  severe  thunderstorms  is  reviewed.  Com- 
prehensive descriptions  of  thunderstorms  and  the  lightning  process  may 
be  found  in  works  written  by  Vi emeis ter , 1  *  **  Malan,^  and  Uman. ^ 

1 . Early  Lightning  Research 

Sir  Isaac  Newton  and  another  English  scientist,  William  Wall, 

were  among  the  first  to  write  of  the  resemblance  between  electricity 

and  lightning.  But  the  electrical  nature  of  thunderstorms  was  not 

verified  until  Benjamin  Franklin's  celebrated  kite  experiment  of  June 
* 

1752.  Besides  the  practical  invention  of  the  lightning  rod,  few 
major  advances  in  the  understanding  of  lightning  were  made  over  the 
next  century.  In  1879  Thomas  A.  Edison  invented  the  light  bulb,  and 
by  1882  the  first  electric  power  lines  had  been  built.  These  were 
often  struck  (or  narrowly  missed)  by  cloud- to-ground  lightning  flashes. 
The  resulting  current  surges  propagated  along  the  power  lines  and 
damaged  equipment.  Power-line  engineers  sought  to  measure  the  peak 
and  duration  of  these  current  surges  in  order  to  devise  methods  of 
neutralizing  their  destructive  effects.  Unfortunately,  equipment  for 
studying  such  transient  phenomena  were  not  then  available. 


* 

D’Alibard  succeeded  in  proving  the  electrical  nature  of  lightning  two 
months  earlier  in  France.  Since  the  experiment  he  performed  had  been 
designed  and  published  by  Franklin,  the  discovery  has  been  attributed 
to  the  latter. 


In  1922  Norinder  adopted  DuFour's  invention  of  the  cathode- 
ray  oscillograph  to  photograph  lightning-induced  electrical  transients. 
During  the  next  decade  this  device  was  used  to  obtain  an  understanding 
of  the  transient  nature  of  lightning.  With  this  knowledge  the  lightning 
arrester  was  perfected  for  power-line  protection. 

The  next  series  of  contributions  can  be  traced  back  to  the 
research  of  Schonland  and  his  associates  in  South  Africa.  In  1933 
they  modified  Boys'  invention  of  the  rotating  camera  to  successfully 
obtain  optical  records  of  lightning  with  extremely  fast  time  resolution. 
Over  the  next  15  years,  analysis  of  their  photographic  records  coupled 
with  simultaneous  electromagnetic  measurements  formed  the  basis  of  our 
current  understanding  of  the  cloud- to-ground  lightning  discharge. 

2 . Previous  Work  on  the  Study  of  Atmospherics 

Even  before  the  invention  of  the  wireless,  Marconi  was  aware 
of  radio  noise  caused  by  lightning;  for  the  Russian  meteorologist  Popov 
had  used  an  invention  remarkably  similar  to  Marconi's  to  receive  atmo- 
spherics for  studying  thunderstorms . The  earliest  research  using  atmo- 
spherics to  understand  thunderstorm  electrification  was  performed  by 
Watson-Watt^ * and  Appleton.1*  Their  studies  were  conducted  at  ELF 
where  atmospherics  are  principally  associated  with  the  return  strokes 
in  a cloud- to-ground  discharge.  These  atmospherics  can  be  characterized 
as  very  strong  pulses  lasting  about  150  ^is. 

Schonland  ° recognized  that  optical  measurements  were  insuf- 
ficient to  characterize  totally  the  electrical  phenomena  of  lightning. 

Accordingly,  he  took  simultaneous  photographic  and  electrical  measure- 

16 

ments  of  close  lightning,  tlis  associate,  Malan,  was  the  first  to 
attempt  to  correlate  the  structure  of  close  atmospherics  at  different 
frequencies.  He  found  that  at  frequencies  above  VLF  the  structure  of 
the  atmospheric  waveform  changed.  The  return  stroke  was  no  longer  the 
dominant  source  of  radiation.  In  fact,  at  HF  a 10-to-15-ms  quiet  period 
followed  each  return  stroke.  During  the  res  of  the  flash,  HF  noise 
radiated  nearly  continuously,  and  intra-cloud  discharges  radiated  as 
much  noise  as  cloud- to-ground  flashes. 


■1 


Horner 


made  a more  comprehensive  study  of  close-by  atmo- 
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spherics  at  different  frequencies.  In  particular  he  showed  that,  in 
disagreement  with  current  theory,  HF  atmospheric  waveforms  could  not  be 
accounted  for  strictly  by  return  strokes  and  stepped  leader  processes. 

He  suggested  that  the  in-cloud  processes  occurring  during  the  intervals 
between  return  strokes  were  the  cause  of  most  HF  radiation. 

19 

Using  receiver  bandwidths  of  over  100  kHz,  Oetzel  and  Pierce 
found  that  HF  atmospherics  were  comprised  of  several  hundred  to  over  a 
thousand  impulses.  The  bandwidth  of  typical  HF  communications  receivers 
is  too  narrow  to  resolve  this  impulse  train.  Thus  at  the  output  of  such 
a receiver  an  HF  atmospheric  usually  takes  the  form  of  a nearly  continu- 
ous burst  of  noise. 

During  the  1950s,  research  was  conducted  to  develop  statisti- 
cal descriptions  of  atmospheric  noise  for  designing  modern  communica- 
tions equipment  and  modulation  formats.  Since  atmospherics  often  pro- 
duce the  dominant  noise  source  below  20  MHz,  HF  noise  was  included  in 
this  study.  The  fact  that  individual  HF  atmospherics  often  overlapped 

and  thus  could  not  be  distinguished  was  of  no  concern  here,  since  the 

20 

object  was  to  characterize  the  noise  background.  Crichlow  et  al . 

developed  an  empirical  model  of  the  amplitude  probability  distribution 

of  atmospheric  noise  received  on  an  omnidirectional  antenna.  Others 

have  used  the  slightly  less  accurate  but  more  tractable  lognormal 

distribution  to  characterize  atmospheric  noise.  Various  statistical 

models  of  the  lightning  radiation  process  were  developed  for  describing 

the  temporal  character  of  atmospheric  noise.  Of  these,  the  impulsive 

21 

noise  model  of  Hall's  has  been  considered  accurate,  particularly  for 

22 

VbF.  The  Poisson-Poisson  model  of  Furutsu  and  Ishida  is  also  appro- 
priate for  modeling  atmospheric  noise  at  HF. 

15 

Schonland  realized  the  potential  benefits  of  receiving 
remote  atmospherics.  From  a single  station,  information  could  be 
gathered  from  many  storms  at  various  distant  locations.  However, 
until  the  present  work,  remote  measurement  programs  have  concentrated 
on  VLF  atmospherics.  Because  of  their  strength  and  their  ability  to 


propagate  over  long  distances  without  great  attenuation,  they  can  be 


easily  received.  And  because  they  take  the  form  of  a few  impulses  per 
flash,  chances  remain  high  that  individual  impulses  of  atmospherics 
from  different  storms  will  not  overlap.  The  principal  source  of  in- 
formation derived  from  these  measurements  has  been  the  location  of 
remote  thunderstorms . Using  a network  of  VLF  direction  finders,  the 
British  have  been  successful  in  monitoring  the  location  of  thunder- 
storms over  portions  of  northern  Europe  and  the  North  Atlantic. 

Horner'*  pointed  out  that  the  structure  of  11F  atmospheric 

noise  would  become  more  variable  when  received  on  a directional  antenna. 
23 

Ortenberger  used  this  principle  to  show  that  the  level  of  atmospheric 

noise  received  on  a directional  antenna  could  vary  substantially  with 

9 

direction.  To  the  author  s knowledge,  Horner  made  the  only  previous 
attempt  to  infer  specific  information  from  individual  HE  atmospherics 
propagated  via  skywave.  Receiving  the  same  atmospheric  at  different 
HF  frequencies  he  was  able  to  use  predictions  of  propagation  conditions 
to  estimate  the  range  to  the  flash  to  an  accuracy  of  100  km. 

3 . Atmospherics  Associated  with  Severe  Thunderstorms 

With  the  appropriate  meteorological  conditions,  thunderstorms 
can  develop  giant  proportions.  These  thunderstorms  are  classified  as 
severe  if  they  produce  at  least  one  of  the  following:  a tornado,  funnel 

cloud,  waterspout,  hail  (greater  than  0.75  inch  diameter),  wind  gusts 
(greater  than  50  knots),  or  extreme  upper  air  turbulence.  The  intensely 
active  electrical  nature  of  severe  storms  has  often  been  noted  and  un- 

2 

usual  lightning  has  even  been  observed  within  the  actual  tornado  vortex. “ 
25 

In  1950,  Jones  reported  that  VLF  atmospherics  radiated  from  severe 

thunderstorms  exhibit  identifiable  characteristics,  thus  suggesting  a 

potential  mechanism  for  severe  thunderstorm  warnings. 

Under  Project.  Tornado-Spherics"^  the  Air  Force  studied  10-kllz 

atmospherics  and  observed  unusually  high  numbers  of  atmospherics  about 

one  hour  prior  to  tornado  sightings.  Typically,  the  rate  of  atmospherics 

27 

returned  to  normal  prior  to  tornado  formation.  Jones  found  that  atmo- 
spheric rates  at  150  kHz  were  a more  sensitive  indicator  of  tornado 
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m 
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activity.  Rates  at  this  frequency  were  unusually  high  before  and  during 
known  tornadoes.  He  also  observed  a localized  region  of  intense  light- 
ning activity  within  the  parent  cloud  and  suggested  that  this  "pulse 

generator"  was  the  origin  of  this  unusual  electrical  activity.  On  the 
28 

other  hand,  Scouten  found  atmospherics  received  on  a 150-kHz  direc- 
tion finder  to  be  distributed  over  a wider  region  within  the  parent 

cloud  centered  about  the  vortex. 

29 

Silberg  studied  the  radiation  spectrum  of  atmospherics  from 

severe  thunderstorms  and  found  it  did  not  decrease  with  frequency  as 

fast  as  for  atmospherics  from  normal  thunderstorms.  In  1968  Weller 

published  a news  article  suggesting  that  home  TV  sets  could  be  used  to 
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warn  of  close  tornadoes.  Biggs  and  Waite  studied  the  cause  of  TV- 

tornado  detections  and  found  a significant  VHF  component  in  atmospherics 
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from  close  tornadoes.  Stanford  and  Lind  further  investigated  severe 
thunderstorm  atmospherics  over  a range  of  frequencies.  They  found  fre- 
quencies above  1 MHz  to  be  most  responsive  to  severe  storm  activity. 

33 

In  the  early  1970s  Taylor  correlated  tornado  activity  with 

the  rate  of  bursts  of  atmospheric  impulses  for  radio  frequencies  ranging 

from  VLF  to  HF.  In  order  to  compare  observations  of  these  different 

frequencies  he  had  to  account  for  the  fact  that  the  number  of  impulses 

radiated  by  a lightning  flash  increased  with  increasing  frequency.  At 

VLF  a normal  cloud-to-ground  discharge  will  produce  about  ten  impulses, 

4 

while  the  same  flash  may  radiate  as  many  as  10  impulses  at  HF.  Using 
fast  logic  circuitry,  Taylor  counted  the  rate  of  impulses  at  different 
amplitude  thresholds  for  frequencies  between  10  kHz  and  3 MHz.  When  an 
impulse  rate  momentarily  exceeded  the  receiving  frequency  measured  in 
kilohertz  a "burst"  was  counted.  Taylor  found  that  at  frequencies 
greater  than  1 MHz  a burst  rate  about  20  per  minute  suggested  the  pres- 
ence of  a severe  thunderstorm . 

34 

Based  on  these  results,  Taylor  constructed  several  atmo- 
spherics receivers  to  record  the  3-MHz  burst  rates  of  storms.  These 
units  were  deployed  throughout  the  Great  Plains  during  the  1972  and 
1973  tornado  seasons.  While  the  correlation  between  tornadoes  and 
high  burst  rates  was  strong,  there  was  also  a disturbingly  high 
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false-alarm  rate.  This  was  partially  attributed  to  the  counting  of 
atmospherics  from  all  storms  within  the  reception  area.  In  an  effort 
to  reduce  this  problem,  some  directional  discrimination  was  provided. 
With  this  new  equipment,  burst  rates  for  eight  45°  azimuth  sectors 
could  be  calculated  independently.  Preliminary  analysis  of  data  re- 
corded using  this  technique  shows  a substantial  reduction  in  false 
alarms  with  no  change  in  detection  rates. 

D . Approach  Used  in  the  Present  Study 

The  use  of  HF  atmospherics  to  study  remote  meteorological  processes 
has  previously  been  neglected.  Propagation  uncertainties  and  confusion 
arising  from  the  reception  of  overlapping  events  combined  to  discourage 
their  use.  In  the  research  described  here,  a directional  HF  antenna  and 
on-line  ionospheric  soundings  were  used  to  diminish  these  difficulties. 

The  initial  phase  of  this  research  was  directed  toward  verifying 
that  convolution  of  the  WARF  antenna  pattern  with  a typical  summertime 
atmospheric  noise  environment  would  result  in  the  strong  reception  of 
remote  atmospherics.  A field  experiment  was  conducted  to  obtain  experi- 
mental verification  that  atmospherics  received  on  the  WARF  antenna  were 
from  lightning  at  a one-hop  distance.  To  the  author's  knowledge,  re- 
sults from  this  test  represent  the  first  time  the  sources  of  remote  HF 
atmospherics  were  positively  located  through  simultaneous  electrostatic 
field  measurements  taken  within  line  of  sight  of  the  lightning. 

Subsequent  emphasis  was  placed  on  the  study  of  HF  atmospherics 
from  remote  severe  thunderstorms.  Data  gathered  during  this  phase  of 
the  research  revealed  a previously  undiscovered  class  of  HF  atmospherics 
associated  with  severe  thunderstorms.  A technique  was  developed  to 
extract  this  component  from  other  types  of  atmospheric  waveforms.  The 
results  were  statistically  analyzed,  confirming  the  correlation  between 
the  reception  of  these  atmospherics  and  the  presence  of  severe  thunder- 
storm events.  These  atmospherics  were  next  compared  with  line-of-sight 
measurements  of  HF  atmospherics  associated  with  severe  storms.  Finally, 
their  temporal  characteristics  were  studied  and  the  results  were  used 
to  postulate  what  meteorological  factors  might  be  their  cause. 
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Original  Contribut ions 


C . 


The  research  tioscribed  here  makes  original  contributions  to  the 
study  of  remote  severe  thunderstorm  activity  and  to  the  understanding 
of  I IF  atmospherics  received  on  directional  antennas.  Specifically: 

(1)  A new  class  of  HF  atmospherics  was  discovered  that  were 
statistically  correlated  with  the  occurrence  of  remote 
severe  thunderstorms.  This  class  of  atmospherics  has  not 
previously  been  isolated  by  line-ol'-sight  measurements 
although  their  short  duration  (lasting  for  less  than  the 
receiver  response  time  of  0.7  ms)  make  them  easily  dis- 
tinguishable front  bursts  of  HF  impulses  associated  with 
atmospherics  from  lightning  in  normal  storms.  Based  on 

a 12-dB  SNK  threshold,  the  rate  of  these  severe  storm 
atmospherics  was  most  commonly  observed  at  120  per  minute 
during  periods  of  severe  storm  activity.  In  contrast,  a 
rate  of  20  per  minute  was  commonly  observed  for  periods 
without  severe  activity.  The  temporal  distribution  of 
these  atmospherics  was  found  to  approximate  that  of  a 
Poisson  process. 

(2)  It  was  shown  that  individual  HF  atmospherics  from  remote 
lightning  can  be  received  on  a directional  antenna  with 
sufficient  strength  and  detail  to  be  used  to  study  these 
remote  lightning  flashes.  For  example,  the  type  of  flash 
(e.g.,  cloud- to- ground  or  intra-cloud)  and  the  rate  of 
occurrence  of  return  strokes  were  measured  from  the 
temporal  properties  of  these  atmospherics. 

Also,  variations  in  the  direction  of  arrival  of  individual 
atmospherics  were  used  to  infer  the  horizontal  dimension 
of  the  region  within  a storm  cloud  that  is  electrically 
active  during  a flash.  For  a storm  occurring  over  New 
Mexico,  the  average  horizontal  dimension  of  the  electrically 
active  region  was  found  to  be  9 + I km. 


II  ANTENNA  FACILITIES  AND  EQUIPMENT 
USED  IN  THE  PRESENT  RESEARCH 


A.  Directional  Antenna 

The  principal  tool  used  throughout  this  research  was  a directional 

HF  antenna — part  of  the  Wide  Aperture  Research  Facility  (WARF).  This 

antenna,  shown  in  Figure  1,  is  a 2.5-km  linear  broadside  array  consisting 

of  256  vertical  monopoles.  As  seen  from  Figure  2,  the  array,  situated 

in  northern  California,  can  be  steered  ±32°,  in  0.25°  steps,  from  its 

boresite  direction  of  true  east  (or  west,  since  it  is  fully  bidirec- 
* 

tional).  Steering  is  achieved  by  electronically  adding  varying  lengths 
of  cable  into  elemental  paths  of  the  analog  beamforming  apparatus.  The 
antenna's  beamwidth  and  directivity  vary  with  frequency  and  steer  angle. 

A nominal  3-dB  beamwidth  of  0.5°  is  achieved  for  a boresite  steer  and 
a midband  frequency  of  15  MHz.  A plot  of  the  theoretical  antenna  pat- 
tern for  the  above  values  is  reproduced  in  Figure  3. 

This  antenna  as  well  as  a close-by  single  vertical  monopole  were 
used  for  recording  HF  noise.  The  data  were  received  on  equipment  built 
at  SRI  for  general  use  in  WARF  experiments.  Data  for  the  experiments 
described  here  were  recorded  on  analog  and  digital  magnetic  tape.  A 
time  code  derived  from  the  site's  cesium-beam  clock  (with  stability  of 
one  part  in  lO*’*  ) was  also  recorded.  Portions  of  the  analog  taped  data 
were  digitized  and  ail  subsequent  processing  was  performed  on  a 16-bit 
Hewlett-Packard  Model  2100  minicomputer. 

Further  specifications  of  this  array  and  the  various  receiving 
and  recording  instruments  available  at  WARF  can  be  found  in  Refs.  3 
and  4.  Since  the  wide-sweep  backscatter  sounding  capability  was  heavily 
used  in  support  of  the  current  research,  it  is  briefly  described  in 
Appendix  A of  this  dissertation. 


* 

Modifications  just  completed  at  the  time  of  this  writing  now  provide 
some  backlobe  suppression. 
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FIGURE  1 THE  WARF  2.5-km  HF  RECEIVING  ANTENNA 
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B. 


Remote  Field  Site 


For  simultaneous  skywave  and  line-of-sight  recordings  of  atmo- 
spherics, a field  site  was  used  in  conjunction  with  the  WARF.  The 
field  site  was  deployed  to  various  parts  of  New  Mexico  to  record  atmo- 
spherics from  close  lightning.  Field-site  equipment  consisted  of  HF 
and  ELF  antennas,  receivers,  and  an  analog  tape  recorder.  A block 
diagram  of  the  equipment  configuration  is  shown  in  Figure  4.  The 
timing  standard  was  synchronized  to  an  accuracy  of  20  ms  by  receiving 
National  Bureau  of  Standards  timing  signals  from  WWV — Boulder,  Colorado. 
Site-to-WARF  synchronization  of  better  than  1 ms  was  later  achieved  by 
adjusting  the  time-base  of  data  received  at  the  remote  site  so  that  the 
first  strong  atmospheric  received  at  both  sites  appeared  coincident. 

The  rough  initial  synchronization  was  sufficient  to  enable  this  atmo- 
spheric association  process  to  work  without  ambiguity.  The  l-part-in-10 
stability  of  the  crystal  standard  used  at  the  remote  site  ensured  that 
the  resulting  adjustment  was  adequate  for  the  remaining  data  collected 
during  that  thunderstorm. 


HF 

ANTENNA 
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FIGURE  4 BLOCK  DIAGRAM  OF  FIELD-SITE  EQUIPMENT  CONFIGURATION 
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The  ELF  receiving  system  was  of  new  construction  and  is  worthy  of 

further  comment.  Its  antenna,  shown  in  Figure  5,  was  reproduced  from 
35 

a design  by  Barham.  The  electronics  consisted  of  a preselector  filter 
and  amplifier  housed  at  the  base  of  the  antenna  and  a receiver  with 
tuned  filters  and  amplifiers. 

The  receiver  was  equipped  with  two  filters  of  different  bandwidths-- 
0.1  to  100  Hz,  and  100  to  1000  Hz.  While  the  first  was  preferable  be- 
cause it  responded  principally  to  changes  in  the  electrostatic  component 
of  close  lightning,  it  could  not  be  used  near  commercial  power  lines. 

At  times  when  such  operations  were  unavoidable  the  100-to-1000-Hz  filter 
was  used. 
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Ill  CHARACTERISTICS  OF  HF  ATMOSPHERIC  NOISE 
RECEIVED  ON  A DIRECTIONAL  ANTENNA 


To  determine  if  a directional  HF  antenna  such  as  that  at  WARF  can 
be  a useful  tool  for  observing  individual  atmospherics  of  remote  origin, 
it  is  necessary  to  investigate  the  interaction  between  antenna  pattern 
and  atmospheric  noise.  At  HF,  non-manmade  noise  is  comprised  of  cosmic 
and  atmospheric  components.  While  the  amplitude  and  time  distributions 
of  cosmic  noise  approach  those  of  thermal,  this  is  not  true  for  atmo- 
spheric noise.  And  such  noise  often  predominates  at  frequencies  below 
20  MHz . 

36 

Beckmann  explained  that  HF  atmospheric  noise  approaches  a 
Rayleigh  distribution  for  small  amplitudes  (and  high  probabilities). 

This  is  presumably  caused  by  the  averaging  or  overlapping  of  indepen- 
dent atmospherics.  For  the  large  amplitudes  (and  low  probabilities) 
the  distribution  diverges  from  Rayleigh,  becoming  more  variable.  These 

large  amplitudes  are  presumably  caused  by  infrequent,  strong  at mo- 

20 

spherics;  perhaps  from  close-by  lightning.  Crichlow  used  a geometric 
construction  of  two  straight  lines  connected  by  a circular  arc  to  des- 
cribe essentially  the  same  cumulative  probability  distribution.  It  is 
significant  that  the  data  used  in  all  these  analyses  were  received  on 
omnidirectional  antennas. 

Horner1  observed  that  the  properties  of  atmospheric  noise  may  be 
dependent  on  the  directivity  of  the  antenna  used  in  their  reception. 

The  noise  becomes  more  variable  with  increasing  directivity — being 
principally  controlled  by  in-beam  thunderstorm  activity.  From  this  we 
might  postulate  that  an  appropriate  model  for  atmospheric  noise  received 
with  a directional  antenna  would  be  composed  of  a low-amplitude  (high- 
probability)  Rayleigh  component  caused  by  the  many  atmospherics  entering 
through  the  antenna's  sidelobes,  and  a comparatively  high- amplitude 
(low-probability)  component  caused  by  in-beam  atmospherics.  As  long 
as  the  storm  is  beyond  the  skip  distance  for  the  frequency  of  observa- 
tion, these  in-beam  atmospherics  need  not  be  close-by  as  for  the  case 
of  an  omnidirectional  antenna.  This  suggests  that  for  a directional 
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antenna  such  as  that  at  WARF,  the  received  noise  will  have,  superimposed 
on  a continuum  of  weaker  Rayleigh  noise,  strong  and  commonly  nonover- 
lapping atmospherics  preserving  enough  detail  to  make  useful  measure- 
ments of  remote  lightning  flashes. 

A . Experimental  Observations  of  Noise  Characteristics  Versus  Antenna 
Directivity 

To  test  the  effect  of  antenna  directivity  on  the  characteristics 
of  received  HF  noise,  data  were  recorded  using  the  0.5°  WARF  antenna 
and  an  omnidirectional  vertical  monopole  situated  about  100  m away.  HF 
noise  was  simultaneously  received  on  two  matched  receivers  of  1.4  kHz 
bandwidth.  The  data  were  recorded  on  a summe”  afternoon  when  no  local 
storm  activity  was  present.  The  WARF  antenna  was  pointed  in  the  general 
direction  of  known  midwestern  thunderstorm  activity. 

Each  receiver  output  was  sampled  at  a 5000-Hz  rate.  The  rms  value, 
E , comPutec*  over  five  samples  (corresponding  to  1 ms)  was  used  to 
estimate  the  noise  envelope.  These  1-ms  samples  were  calculated  for  a 
30-s  interval  of  data  and  used  to  construct  empirical  amplitude-proba- 
bility distributions.  Representative  samples  of  the  resulting  distribu- 
tions are  shown  in  Figure  6.  These  cumulative  distributions  have  been 
plotted  on  Rayleigh  probability  paper,  which  has  the  characteristic  that 
a Rayleigh  distribution  takes  the  form  of  a straight  line  of  slope  -0.5. 
The  distributions  of  Figure  6 exhibit  low- amplitude  Rayleigh  components 

with  more  variable  high- amplitude  components  in  agreement  with  the 

36  20 

distributions  observed  by  Beckmann  and  Crichlow. 

The  distributions  in  Figure  6(a)  have  been  normalized  by  the  rms 

average  of  their  envelopes  in  accordance  with  the  conventions  used  by 

the  authors  mentioned  above.  However,  for  the  purpose  of  comparing  the 

distributions  obtained  from  different  antennas,  it  is  convenient  to 

equalize  the  common  Rayleigh  noise  "background."  This  was  accomplished 

in  Figure  6(b)  by  renormalizing  the  distributions  by  their  respective 
* 

mode  values. 

* 

The  mode  represents  the  value  for  which  the  amplitude-probability 
density  function  attains  its  maximum. 


j 
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PERCENTAGE  OF  TIME  ORDINATE  IS  EXCEEDED 
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FIGURE  6 CUMULATIVE  AMPLITUDE-PROBABILITY  DISTRIBUTIONS  OF  NOISE  RECEIVED 
ON  ANTENNAS  OF  DIFFERENT  DIRECTIVITY 


In  Figure  7 the  data  of  Figure  6 have  been  replotted  as  amplitude- 
probability  density  functions  on  a linear  scale.  This  figure  demon- 
strates that  the  mode  is  strictly  a function  of  the  Rayleigh  portion 
of  the  distribution.  Other  statistics  such  as  the  rms,  mean,  and  to 
a lesser  extent,  median,  are  dependent  on  the  variable  high- amplitude 
component  of  these  distributions. 

From  Figure  6(b)  it  is  evident  that  the  high- amp litude  component 
of  the  observed  atmospheric  noise  is  more  pronounced  for  the  directional 
antenna.  Review  of  the  data  in  the  time  domain  is  necessary  to  deter- 
mine if  this  component  is  produced  by  individual  atmospherics.  Figure  H 
plots  one-second  intervals  of  the  data  used  to  compute  the  distributions 
of  Figure  6.  The  envelope  of  the  noise  received  with  the  two  different 
antennas  has  been  normalized  by  their  respective  mode  values  to  provide 
equal  noise  backgrounds . It  can  be  seen  that,  most  of  the  hi  gh-amp  litude 
values  are  clustered  in  discrete  bursts.  Although  hardly  observable  on 
the  monopole  data,  the  time/amplitude  characteristics  of  the  noise  bursts 
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received  on  the  directional  antenna  are  similar  to  observations  of 

16  18 

single  HF  atmospherics  radiated  from  close  lightning  flashes.  ’ 

(The  receiver  passband  was  too  narrow  to  resolve  individual  impulses 

and  the  atmospheric  therefore  takes  the  form  of  a continuous  burst  of 
22 

noise.)  Furutsu  showed  that  HF  noise  distributions  with  such  un- 
usually strong  amplitude  components  can  be  caused  when  local  thunder- 
storms are  present.  But,  since  no  local  storms  were  present,  the  weight 
of  evidence  suggests  that  these  strong  atmospherics  originated  from 
remote  flashes  (presumably  at  a one-hop  distance)  and  were  received  via 
skywave  propagation. 

An  indirect  method  of  testing  this  conclusion  was  to  estimate  a 
likely  thunderstorm  spatial  distribution  for  one-hop  distances  from 
the  WARF  antenna  and  determine  if  the  resulting  atmospheric  noise  is 
comparable  to  the  noise  observed. 

23 
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A Model  for  One*- Hop  Atmospheric  Noise  Keep  1 ved  at  WARF 


To  construct  a model  of  atmospheric  noise  incident  on  the  WARF 
antenna,  it  is  first  necessary  to  describe  the  thunderstorm  activity 
likely  to  create  that  noise. 

The  flashing  rate  of  a thunderstorm  varies  over  the  course  of  eacli 
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storm  and  from  storm  to  storm.  Cianos  notes  that  observations  of 

flashing  rates  have  varied  from  one  to  more  than  ten  flashes  per  minute 

with  a mean  of  three.  The  duration  of  each  flash  also  varies  with  most 

37 

measurements,  showing  an  average  of  500  ms.  Tfie  atmospheric  noise 

23 

received  on  an  antenna  can  be  modeled,  after  Ortenburger,  as  the 
summation  of  atmospherics  from  many  storms--each  contributing  an  amount 
weighted  by  its  propagation  path  and  the  antenna  pattern. 

Near  the  WARF  receiving  antenna  located  in  the  Central  Valley  of 
northern  California,  there  are  seldom  any  local  thunderstorms.  The 
Sierra  Nevada  Mountain  Range  provides  a moderate  source  of  line-of-sight 
storms,  but  groundwave  losses  over  this  80-km  path  are  comparable  to 
one-hop  skywave  losses.  As  a consequence,  the  major  source  of  HF  atmo- 
spheric noise  received  at  WARF  propagates  via  ionospheric  reflection. 
Because  of  lower  path  losses,  those  storms  in  the  midwest  having  a one- 
hop  propagation  path  will  predominate  over  storms  from  the  more  elec- 
trically active  tropical  regions,  which  require  multiple-hop  propagation 
paths . 

Figure  9 shows  a representative  example  of  the  region  having  strong 

one-hop  propagation  to  WARF.  This  region  will  change  with  choice  of 

receiving  frequency  and  ionospheric  conditions;  however,  it  is  useful 

for  estimating  the  approximate  azimuthal  f lashing- rate  density.  To 

simplify  calculations  we  will  assume  propagation  losses  to  be  constant 

* 

throughout  the  one-hop  region.  The  numbers  shown  in  Figure  9 are 


* 

The  number  of  ocean  thunderstorms  is  thought  to  be  negligible  compared 
to  land  storms.  Recent  WARF  antenna  modifications  provide  east-west 
backlobe  discrimination,  and  initial  measurements  support  the  assumption 
that  most  atmospheric  noise  received  at  WARF  originates  to  the  east  of 
the  antenna. 
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FIGURE  9 AVERAGE  MONTHLY  SUMMERTIME  THUNDERSTORM  — DAYS  WITHIN  A 
TYPICAL  ONE-HOP  PROPAGATION  REGION  FROM  WARF 

T 8 

derived  from  the  summertime  monthly  thunderstorm-day  statistic/  These 
statistics  estimate  the  number  of  days  per  month  that  thunder  is  heard 
at  a particular  location. 

S9 

Using  a relationship  developed  by  Pierce,  we  can  estimate  the 

flashing-rate  density  from  the  thunderstorm-day  statistic.  If  we 

assume  an  average  of  9 for  the  summertime  monthly  thunderstorm-day 

statistic  in  the  one-hop  region  of  Figure  9,  the  afternoon  flashing- 
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rate  density  can  be  estimated  at  1.8  x 10  km  -min.  Applying  this 
spatial  flashing-rate  density  to  the  geometry  of  Figure  9 results  in 
an  average  azimuthal  flashing-rate  density  of  1.1  flashes  per  degree- 
minute.  Over  the  electrically  active  one-hop-by-60°  sector  shown  in 
Figure  9 (and  ignoring  the  remaining  relatively  nonactive  one-hop 
region),  this  average  flashing  rate  would  yield  a total  of  265  flashes 


* 

per  minute.  Since  each  flash  lasts  several  hundred  milliseconds , this 
calculation  indicates  that  radiation  from  several  one-hop  flashes  may 
be  received  simultaneously.  Thus,  for  an  omnidirectional  antenna  the 
received  atmospheric  noise  will  be  nearly  continuous  and  atmospherics 
from  individual  flashes  will  not  be  easily  recognized.  In  contrast, 
for  the  0.5°  antenna  at  WARF,  the  average  in-beam  flashing  rate  is  only 
2.2  flashes  per  minute.  Because  this  narrow  beam  intersects  a rela- 
tively small  area,  localized  weather  conditions  make  it  likely  that  the 
actual  in-beam  flashing  rate  will  substantially  deviate  from  the  average. 
However,  even  considering  an  in-beam  flashing  rate  as  high  as  10  per 
minute,  most  flashes  will  be  separated  from  one  another  in  time.  If 
the  antenna  gain  and  sidelobe  levels  are  adequate,  these  in-beam  events 
should  be  strong,  compared  to  the  continuum  noise  background  produced 
by  the  many  out-of-beam  events. 

In  conclusion,  it  should  be  stated  that  the  above  model  of  HF 
atmospheric  noise  incident  at  WARF  must  be  considered  as  a very  rough 
approximation.  Propagation  or  localized  weather  conditions  may  easily 
place  this  model  in  error  by  an  order  of  magnitude.  However,  the  cal- 
culations are  useful  in  substantiating  the  preliminary  conclusion  that 
the  directional  antenna  at  WARF  has  the  capability  to  observe  atmo- 
spherics from  individual  remote  lightning  flashes. 


* 

It  is  expected  that  these  flashes  would  not  be  evenly  distributed  over 
such  a large  area.  However,  only  the  total  flashing  rate,  not  its 
azimuthal  distribution,  is  of  importance  when  calculating  the  atmo- 
spheric noise  received  by  an  omnidirectional  antenna  or  through  the 
sidelobes  of  a directional  antenna. 
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IV  OBSERVATIONS  OF  REMOTE  HF  ATMOSPHERICS 


A.  Verification  of  the  Source  of  Atmospherics  Observed  with  the 

Directional  Antenna 

The  observations  described  in  the  previous  section  imply  that  many 
of  the  strong,  and  well  defined  atmospherics  received  on  the  WARF  an- 
tenna originate  at  a remote  distance  and  then  propagate  to  the  antenna 
via  skywave.  Direct  verification  of  this  hypothesis  was  made  by  field 
measurements.  An  experiment  was  designed  to  record  and  compare  signals 
received  on  the  0.5°  antenna  with  measurements  simultaneously  recorded 
in  the  immediate  vicinity  of  a thunderstorm  at  a one-hop  distance  from 
northern  California. 

1 . Design  of  Experiment 

A remote  field  site  was  established  in  an  area  of  New  Mexico 
likely  to  produce  thunderstorm  activity.  When  a storm  approached  this 
site,  HF  noise  data  were  simul taneously  recorded  there  and  in  California 
with  the  0.5°  antenna  pointed  in  the  storm's  direction.  Wide-sweep 
backscatter  soundings  were  taken  at  WARF  (see  Appendix  A)  to  determine 
the  best  band  of  frequencies  available  for  propagation  to  the  remote 
site  and  thereby  infer,  by  reciprocity,  the  best,  frequency  band  avail- 
able from  the  remote  site  to  WARF.  The  1.4-kllz  receivers  of  both  sites 
were  then  tuned  to  an  interference-free  channel  within  this  frequency 
band.  This  selection  process  was  principally  accomplished  by  listening 
to  the  receiver  output  for  the  presence  of  man-made  signals.  Simul- 
taneous HF  records  were  then  collected  at  both  sites  for  use  in  making 
a comparison  of  atmospheric  noise  received  via  direct  and  one-hop  paths. 

ELF  records  were  also  recorded  at  the  field  site.  These 
records  were  intended  for  use  in  identifying  flashes  close  to  the  site. 
The  characteristics  of  ELF  radiation  from  lightning  change  rapidly  with 
propagation  distance.*0  At  close  ranges  the  electrostatic  field  is 
strong.  But  this  field  decreases  in  proportion  to  the  cube  of  distance. 
The  weaker  radiation  field,  which  decreases  linearly  with  distance,  will 
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predominate  at  longer  ranges.  The  reception  of  ELF  radiation  indicative 
of  a strong  electrostatic  component  can  thus  be  used  to  identify  flashes 
occurring  close  to  the  field  site.  The  simple  form  of  the  electrostatic 


field  can  also  be  used  to  identify  specific  events  within  a lightning 
flash.40’11 


2 . The  Remote  Field  Site 

An  area  in  north  central  New  Mexico  was  selected  for  the  loca- 
tion of  the  field  site.  This  region  was  chosen  because  of  its  vigorous 
thunderstorm  activity  and  its  good  one-hop  propagation  path  back  to 
California.  Because  the  ELF  measurement  is  sensitive  to  60-Hz  powerline 
radiation,  actual  field  sites  were  located  several  miles  from  commercial 
power.  The  equipment  was  housed  in  a van  and  powered  by  a gas-driven 
electric  generator.  Two  New  Mexico  sites  were  selected.  The  first  was 
located  near  Bernalillo  and  the  second  was  situated  near  Redondo  Peak 

V 

in  the  Jemez  Mountain  Range. 

This  experiment  was  conducted  on  eleven  afternoons  during  the 
summer  of  1974.  Due  to  unusually  dry  conditions  prevalent  in  New  Mexico, 
only  one  thunderstorm — at  the  start  of  the  first  day  of  operations — 
came  close  to  either  field  site.  During  equipment  setup  on  the  after- 
noon of  31  July  1974  an  unusually  active  thunderstorm  developed  and 
slowly  passed  near  the  Bernalillo  field  site.  Figure  10  shows  a New 
Mexico  radar  precipitation  map  taken  at  2230  UT , 31  July  1974.  The 
field  site  location,  and  3-dB  and  10-dB  footprints  of  the  WARF  antenna 
have  been  added.  While  the  HF  equipment  at  the  field  site  was  not  yet 
operational,  good  ELF  records  were  obtained.  HF  noise  received  with 
the  WARF  antenna  was  simultaneously  recorded  in  California.  These  data 
were  used  to  compare  noise  measurements. 

3 . Data  Comparison 

The  analog  tape  recordings  from  the  New  Mexico  field  site 
were  brought  back  to  California  and  digitized  for  comparison  with  the 
WARF  HF  records.  Figures  11(a)  and  (c)  reproduce  ELF  measurements 


FIGURE  10  FOOTPRINT  OF  DIRECTIONAL  ANTENNA  OVERLAID  ON  WEATHER  RADAR  PRECIPITATION  MAP 


recorded  in  New  Mexico.  Figures  11(b)  and  (d),  respectively,  plot  the 
magnitude  of  HF  noise  simultaneously  received  at  WARF.  The  time  axes 
have  been  shifted  11  ms  to  account  for  differential  propagation  delays 
and  time  synchronization  errors  as  outlined  in  Section  II,  This  shift 
has  the  effect  of  arbitrarily  forcing  a correlation  between  one  feature 
of  the  first  strong  atmospheric  received  at  WARF  and  an  equivalent  fea- 
ture received  at  the  New  Mexico  site.  Although  the  11-ms  shift  biased 
the  probability  of  making  this  first  coi-relation,  results  of  an  analysis 
described  in  Appendix  B proved  that  subsequent  events  were  free  of  bias. 
As  shown  in  Appendix  B,  it  is  extremely  unlikely  that  the  observed  cor- 
relation between  features  of  atmospherics  received  at  both  sites  could 
have  occurred  merely  by  chance. 

Analysis  of  the  ELF  records  indicate  that  they  were  produced 
by  close  lightning  flashes.  The  abrupt  positive  field  change  is  indica- 
tive of  cloud- to-ground  return  strokes,  while  the  slight  positive 

deflection  just  prior  to  the  first  stroke  indicates  the  initial  break- 
down and  stepped  leader  phases.  A comparison  of  these  ELF  records  with 
the  HF  noise  bursts  received  at  WARF  shows  a high  degree  of  correlation. 
The  HF  noise  starts  with  the  onset  of  the  breakdown  process  and  increases 
in  amplitude  until  the  return  stroke  when  it  abruptly  stops  for  a short 
period.  This  quenching  of  HF  noise  has  been  observed  in  records  taken 
close  to  cloud-to-ground  lightning . ^ It  is  characteristic  of  the 

return  stroke  and  can  be  used  for  identification  of  cloud- to-ground 
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discharges . 

Over  a 5-minute  period,  87%  of  the  flashes  recorded  in  New 
Mexico  were  correlated  with  HF  noise  received  on  the  WARF  antenna  in 
California.  For  most  of  these  events,  the  atmospherics  received  at 
WARF  were  sufficiently  strong  to  identify  the  type  of  flash  and  the 
duration  of  the  initial  breakdown  process  for  cloud- to- ground  flashes. 
Radiation  prior  to  the  first  return  stroke  of  a multiple-stroke  cloud- 
to-ground  flash  typically  produced  the  strongest  signal.  Radiation 
from  higher-order  strokes  was  not  always  observed. 
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These  observations  can  be  made  because  the  effects  of  iono- 


spheric propagation  do  little  to  alter  the  shape  of  the  received  atmo- 
spheric. Ionospheric  multipath  and  dispersion  will  change  the  shape 
and  number  of  impulses  contained  in  an  atmospheric.  However,  the 
1.4-kHz  receiver  bandwidth  used  in  this  study  responds  to  the  envelope 
of  this  impulse  train,  which  is  insensitive  to  propagation  effects. 
Other  ionospheric  effects  such  as  amplitude  fading  are  generally  too 
slow  to  alter  individual  atmospherics,  which  typically  last  less  than 
one  second. 

Not  all  atmospherics  received  at  WARF  originated  close  to 
the  New  Mexico  field  site.  Figure  12  plots  one  minute  of  received 
noise  on  a compressed  time  scale.  Periods  when  close  flashes  were  re- 
corded on  the  New  Mexico  equipment  are  underlined.  The  cross-hatched 
period  indicates  atmospherics  received  at  WARF  and  correlated  with  New 
Mexico  records  but  also  showing  evidence  of  other-event  overlap.  For 
the  5-minute  interval  analyzed,  the  Bernalillo  storms  accounted  for  35% 
of  the  atmospherics  received  in  California.  Thunderstorm  activity  over 
the  California  Sierra  Nevada  mountains  and  in  eastern  New  Mexico  (see 
Figure  10)  most  likely  caused  the  remaining  atmospherics. 

While  it  was  beyond  the  scope  of  this  experiment  to  account 
for  all  atmospherics,  the  data  directly  verified  that  strong  atmo- 
spherics of  remote  origin  can  be  received  on  a directional  antenna. 
Furthermore,  event  overlapping  is  sufficiently  reduced  that  individual 
atmospherics  can  be  observed  and  used  to  identify  events  within  a flash. 

B.  Variations  in  the  Bearing  of  Noise  Received  from  a Single 

Lightning  Flash 

The  previous  subsection  established  that  information  contained  in 
the  amplitude-versus-time  waveform  of  atmospherics  can  be  remotely 
monitored  using  a directional  antenna.  In  this  subsection  we  will  show 
that  the  atmospheric  wavefront  contains  information  on  horizontal  varia- 
tions of  the  radiation  source  involved  within  a single  remote  flash. 

This  surprising  observation  was  noted  when  the  WARF  was  configured  to 
obtain  accurate  bearing  estimates  of  remote  atmospherics. 
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FIGURE  12  ATMOSPHERIC  NOISE  RECEIVED  WITH  THE  WARF  ANTENNA  POINTED 
IN  THE  DIRECTION  OF  THE  BERNALILLO  STORMS 


1 . Bearing  Estimation  of  Atmospherics 

A digital  beamforming  technique  was  used  with  the  WARF  antenna 
to  produce  eight  simultaneous  0.5°  beams  spaced  0.5°  apart.  This  tech- 

5 

nique  is  similar  to  that  developed  by  Sweeney  but  is  simplified  by  the 
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addition  of  an  eight-channel  coherent  receiver  with  matched  passbands 
and  eight  simultaneous  analog- to-digital  converters.  The  2.5-km  antenna 
is  divided  into  eight  equal  parts.  An  analog  beam  is  formed  for  each 
of  the  eight  subarrays.  Noise  from  each  of  these  4°  beams  is  received 
on  one  of  the  eight  receiver  channels.  The  outputs  are  digitized  and 
processed  through  a two-dimensional  Fourier  transform  to  obtain  bearing 
estimates  of  the  received  noise.  Independent  bearing  measurements  are 
computed  for  every  4 ms  of  data.  (Further  details  of  this  digital 
beamforming  technique  appear  in  Appendix  C.) 
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When  eight  contiguous  beams  are  lormect  simultaneously,  the 
received  noise  energy  versus  azimuth  and  time  can  be  computed.  While 
azimuthal  resolution  is  fixed  by  antenna  length,  single-source  azi- 
muthal accuracy  is  limited  only  by  SNR.  For  skywave  propagation, 
absolute  accuracy  is  further  limited  by  coning  errors  and  by  deviations 
from  great-circle  propagation  due  to  ionospheric  tilts.  These  factors 
do  not  affect  relative  azimuthal  accuracy  (i.e.,  measurements  of  azi- 
muth difference).  With  adequate  SNR,  interpolation  can  be  used  on  the 
digital  beamformed  output  to  measure  the  relative  azimuth-versus-time 
structure  of  an  atmospheric  to  an  accuracy  of  0.06°. 

The  beamformer  output  can  be  described  with  the  aid  of 
Figure  13.  The  noise  received  on  one  of  the  eight  subarrays  is  plotted 
in  Figure  13(a)  in  magnitude-versus-time  format.  Figure  13(b)  plots 
these  data  as  beamformed  noise-magnitude  versus  time  and  azimuth  in  a 
pseudo-three-dimensional  format.  The  atmospheric  shown  in  Figure  13 
was  recorded  on  10  April  1975.  Its  origin  is  unknown  but  was  probably 
a storm  system  over  eastern  New  Mexico  and  the  Texas  panhandle.  The 
azimuthal  variation  of  this  atmospheric  was  confined  to  within  a 0.15° 
sector . 

2 . Observations  of  Azimuthal  Variations 

Figure  14  plots  the  amplitude  versus  time  and  azimuth  of  a 
record  containing  an  atmospheric  recorded  roughly  20  minutes  after  that 
of  Figure  13.  During  the  700  ms  of  this  burst,  the  apparent  direction 
of  arrival  changes  by  0.7°.  Analysis  of  ten  atmospherics  received  on 
10  April  shows  an  average  azimuthal  variation  of  0.5°  ±0.2°. 

A horizontal  variation  in  the  electrically  active  cloud  region 
during  a lightning  flash  is  one  plausible  explanation  of  these  observa- 
tions. Work  by  Ogowa,  **  Proctor,  '->  and  Few  ^ have  recently  shown  that 
the  in-cloud  charge  region  involved  during  a flash  commonly  has  a 
larger  horizontal  component  than  was  previously  assumed.  Because  the 
range  to  the  origin  of  the  10  April  atmospherics  was  unknown,  the 
cross-range  distance  corresponding  to  the  observed  azimuthal  change 
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(a)  SUBARRAY  OUTPUT 


(b)  BEAMFORMED  OUTPUT 


* 

can  only  be  estimated.  Without  direct  verification,  the  possibility 
of  equipment  malfunction  or  ionospheric  propagation  peculiarities  must 
be  considered  before  presuming  that  these  observations  were  actually 
the  result  of  changes  within  the  flash.  A second  New  Mexico  field  test 
ruled  out  the  possibility  of  these  alternative  explanations. 

During  July  of  1975  a field  site  in  Los  Lunas,  New  Mexico  was 
manned  for  an  unrelated  experiment.  The  HF  and  ELF  receiving  equipment 
used  in  the  field  experiment  described  in  the  previous  section  was  in- 
stalled and  at  the  ready  if  a thunderstorm  should  occur  close  by.  Be- 
cause of  the  proximity  to  power  lines,  the  100-to-1000-Hz  ELF  passband 
was  used.  On  10  July  1975  a thunderstorm  passed  close  to  the  Los  Lunas 
site.  Data  from  this  storm  were  simultaneously  recorded  in  New  Mexico 
and  at  WARF. 

Changes  in  the  received  azimuth  of  a single  atmospheric  were 
again  observed  on  10  July.  Figure  15  shows  a one-second  record  of  the 
Los  Lunas  ELF  (a)  and  HF  (b)  channels  along  with  the  subarray  (c)  and 
beamformed  (d)  data  received  simultaneously  at  WARF.  The  horizontal 
line  through  Figure  15(d)  indicates  the  observed  bearing  of  a calibra- 
tion signal  received  at  WARF  and  transmitted  from  Los  Lunas.  During 
the  one  second  of  data  displayed  in  Figure  15,  noise  from  two  flashes 
was  received  with  the  WARF  antenna.  The  first  flash  evidently  origi- 
nated from  a storm  other  than  the  one  over  Los  Lunas.  This  is  concluded 
because  noise  from  this  flash  was  not  received  in  Los  Lunas  and  the 
observed  azimuth  was  0.84°  away  from  the  Los  Lunas  bearing.  In  con- 
trast, the  second  flash  appears  to  have  come  from  the  Los  Lunas  storm. 

The  received  noise  from  this  flash  shows  a total  bearing  change  of 

o o 

0.52  over  120  ms  for  an  average  cross-range  rate  of  4.3  /s. 

* 

Based  on  the  propagation  conditions  observed  on  the  wide-sweep  back- 
scatter  soundings,  atmospherics  from  storms  between  950  and  1600  km 
might  have  been  received.  Weather  Bureau  data  show  only  one  storm 
system  within  this  area.  The  range  to  this  storm  was  approximately 
1500  km,  suggesting  an  average  cross-range  variation  of  13  km. 
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During  data  collection  from  this  storm  a signal  was  periodi- 
cally transmitted  from  the  Los  Lunas  site  and  received  at  WARF.  The 
amplitude  versus  azimuth  and  time  of  this  signal  was  computed  by  the 

same  beamforming  program.  A 0.24°  azimuthal  variation  of  this  refer- 

o 

ence  signal  was  observed  over  27  s with  a maximum  rate  of  0.04  /s. 

This  is  two  orders  of  magnitude  smaller  than  the  rate  observed  in  the 
atmospheric  of  Figure  15. 

The  reference  signal  propagated  via  essentially  the  same 
ionospheric  path  and  was  received  and  processed  in  the  same  way  as 
atmospherics  received  from  the  storm.  Since  azimuthal  fluctuations  in 
the  calibrate  signal  are  hardly  detectable  on  a one-second  time  scale, 
it  can  be  assumed  that  neither  equipment  nor  ionospheric  effects  can 
explain  the  apparent  directional  changes  observed  in  some  atmospherics. 

The  azimuthal  extent  of  33  flashes  determined  to  have  origi- 
nated in  the  Los  Lunas  storm  were  analyzed.  Using  1340  km  as  the  slant- 
range  distance  between  the  WARF  antenna  and  Los  Lunas  (appropriate  for 
the  propagation  conditions  present),  the  cross  ranges  corresponding 
to  the  measured  azimuthal  changes  are  plotted  in  Figure  16.  On  the 

4 

average,  these  values  are  larger  than  the  3 to  8 km  measured  by  Proctor 
(at  VHF)  and  Ogawa 1 * (at  ELF).  On  the  other  hand,  the  distribution  of 
lengths  is  remarkably  similar  to  those  measured  by  Few ^ (using  acoustic 
signals)  for  frontal  type  thunderstorms.  It  should  be  emphasized  that 
the  various  results  cited  here  were  obtained  using  different  measurement 
techniques  and  hence  are  sensitive  to  different  processes  within  the 
flash . 

In  summary,  observations  of  atmospherics  described  here  show 
substantial  changes  in  direction  of  arrival,  which  were  not  observed 
on  a ground-based  calibration  source.  While  the  elevated  geometry  of 
the  lightning  radiation  might  account  for  minor  differences,  it  is  not 
a likely  explanation  for  the  substantial  differences  observed  between 
the  azimuthal  characteristics  of  these  two  signal  sources.  Although 
it  has  not  been  confirmed  by  direct  measurements,  the  data  provide 
strong  evidence  that  the  measured  direction  of  arrival  of  remote  HF 
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V 1 WE  ST  I GAT  I ON  OF  HF  ATMOSPHERICS 

FROM  SEVERE  i ..UNDERSTORMS 


Experiments  described  in  the  previous  sections  demonstrated  the 
ability  of  a directional  HF  antenna  to  observe  temporal  and  spatial 
properties  of  individual  remote  lightning  flashes  through  the  recep- 
tion of  the  atmospherics  they  radiate.  This  capability  was  applied  to 
observations  of  atmospherics  from  remote  severe  thunderstorms.  Over 
the  past  two  decades  researchers  have  noted  that  severe  storms  are 
usually  accompanied  by  increased  electrical  activity.  Increased  radia- 
tion of  atmospherics  during  severe  thunderstorms  has  been  observed  at 

frequencies  between  VLF“'’’“^  and  UHF,'^  with  recent  work  concentrating 
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in  the  HF  band.'  ’ However,  the  author  is  unaware  of  any  attempt, 
prior  to  the  work  described  here,  to  use  atmospherics  to  study  severe 
storms  from  ranges  beyond  line  of  sight. 

In  the  case  of  line-of-sight  measurements,  where  the  spreading 
loss  of  the  radiation  field  predominates,  amplitude  threshold  levels 
can  be  set  so  that  only  close-by  HF  atmospherics  will  be  counted.  For 
this  reason  the  total  count  of  atmospherics  can  be  used  for  sensing 
the  presence  of  local  severe  thunderstorms.  However,  for  remote  mea- 
surements it  is  possible  that  several  thunderstorms  not  associated  with 
the  severe  storm  may  also  lie  within  the  antenna's  main  beam.  In  this 
case  the  total  count  of  received  atmospherics  would  not  be  a sensitive 
measure  of  the  severity  of  any  one  storm.  Thus  for  the  remote  study 
of  severe  thunderstorms  to  be  possible,  the  storms  must  radiate  atmo- 
spherics that  can  be  uniquely  identified  despite  the  presence  of 
extraneous  atmospherics.  This  section  describes  an  experiment  con- 
ducted to  determine  if  such  a class  of  identifiable  severe  storm  atmo- 
spherics exist.  Initial  positive  indications  are  presented.  The  next 

section  will  describe  a statistical  analysis  of  t h>  observed  effect. 

> 
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A . Experimental  Configuration 

Clearly,  for  the  study  of  remote  severe  thunderstorms  via  1H' 
atmospherics,  data  must  be  collected  with  an  antenna  pointed  in  the 


direction  of  such  a storm.  But  the  precise  time  and  location  of  severe 

thunderstorms  cannot  be  predicted.  Their  occurrence  is  highly  sporadic 

and  widely  dispersed  geographically — even  during  the  springtime  peak  ol 

the  tornado  season.  Furthermore,  severe  storms  are  seen  only  through 

their  more  observable  manifestations  (e.g.,  tornadoes,  large  hail, 

etc.),  and  the  brevity  of  these  events  (the  average  lifetime  of  a 
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tornado  is  15  minutes)  coupled  with  the  delays  inherent  in  current 
warning  procedures  make  it  impractical  to  gather  useful  data  using 
reports  of  sightings. 

An  alternative  to  pointing  the  antenna  in  the  direction  of  a known 
severe  storm  is  to  scan  the  antenna  over  a series  of  bearings  likely  to 
contain  such  activity.  In  this  way  the  noise  originating  within  a 
sizable  region  can  be  recorded.  Later,  the  data  can  be  analyzed  to 
discover  possible  effects  caused  by  storms  with  severe  events  sighted 
in  that  region.  As  with  any  scanning  antenna,  the  larger  the  azimuth 
sector  scanned,  the  less  time  noise  is  received  from  any  one  direction. 
If  a reasonable  duty  factor  is  to  be  maintained,  the  size  of  the  sector 
monitored  must  be  limited.  This  was  accomplished  by  using  National 
Weather  Service  severe-weather  forecasts^**’ ^ to  select  those  regions 
within  the  country  most  likely  to  contain  severe  storm  activity. 

The  National  Severe  Storms  Forecast  Center  (NSSFC)  forecasts  areas 
likely  to  develop  severe  thunderstorms.  They  issue  a "Severe  Storm 
Watch"  to  identify  a specific  region  (nominally  110  miles  by  200  miles) 
and  time  period  during  which  tornadoes  are  likely  to  occur.  The  loca- 
tion of  these  regions  may  vary  widely  from  day  to  day  or  even  during 
the  course  of  one  day.  By  the  use  of  these  watches  the  antenna  scan 
sector  can  be  limited  while  a reasonable  probability  is  maintained  of 
receiving  noise  from  the  direction  of  a severe  thunders to rm . 

The  antenna  scanning  formats  chosen  for  this  experiment  represent 
a compromise  among  revisit  time  (the  time  between  successive  looks  in 
the  same  direction),  dwell  time  (the  time  during  which  data  are  re- 
ceived from  one  direction),  azimuthal  step  size  (bearing  difference 
between  successive  beam  positions),  and  azimuth  coverage  (total  azimuth 
sector  scanned).  A five-minute  revisit  time  was  deemed  the  maximum 


allowable,  considering  the  brevity  of  some  severe  storm  events.  Thirty 
seconds  was  considered  the  minimum  dwell  time  for  evaluation  of  the 


received  noise  structure.  Since  most  severe  storm  watches  occupy  an 
azimuth  sector  from  10°  to  20° , the  above  constraints  lead  to  a 2° 
choice  of  antenna  step  size.  Table  1 describes  the  resulting  two 
scanning  formats  used  for  this  experiment.  Selection  between  these 
two  formats  depended  on  current  watch  information. 

Table  1 

ANTENNA  SCAN  FORMATS 


Format  A 

Format  B 

Revisit 

time 

5 min 

5 min 

Dwell  time 

60  s 

80  s 

Azimuth 

step  size 

2° 

2° 

Azimuth 

coverage 

8° 

o 
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Severe  storm  watches  were  communicated  to  the  WARF  from  the  Naval 
Weather  Station  at  Moffett  Field  and  from  the  San  Francisco  branch  of 
the  National  Weather  Service.  At  WARF  the  watch  area  was  plotted  on  a 
U.S.  map  overlayed  with  a grid  of  antenna  bearings  and  approximate  one- 
hop  F-layer  propagation  time  delays.  An  example  is  reproduced  in 
Figure  17.  The  watch  area  is  principally  contained  within  a 10°  sector, 
and  thus  Format  A of  Table  1 was  selected  with  bearing  limits  of  100° 
and  108  . With  the  use  of  wide-sweep  backscatter  soundings,  a radio 
frequency  appropriate  for  receiving  noise  from  the  watch  area  was 
selected  by  the  technique  described  in  Appendix  A.  The  selected  fre- 
quency was  changed  from  time  to  time  as  ionospheric  propagation  and 
man-made  interference  conditions  changed. 

B.  Data  Recording  and  Monitoring 

As  the  antenna  stepped  across  the  scan  sector,  the  noise  output 
from  the  1.4-kHz  receiver  was  analog  tape  recorded  along  with  the 
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FIGURE  17  A.  U.S.  MAP  SHOWING  A SEVERE  STORM  WATCH  AREA  AND  A WARF 
DIRECTION/TIME-DELAY  OVERLAY 


wide-sweep  backscatter  soundings  and  a time  code.  Data  were  collected 
in  this  fashion  on  several  afternoons  in  the  latter  half  of  May  1975. 
During  data  collection  the  noise  channel  was  connected  to  a speaker 
and  continually  monitored  by  ear.  The  frequency  spectrum  of  the  re- 
ceiver output  was  also  plotted  to  provide  a continuous  hard  copy  of 
the  received  data.  These  outputs  were  convenient  for  monitoring  over- 
all atmospheric  noise  activity  as  well  as  weak  interference. 

This  real-time  spectral  record  was  obtained  by  passing  the  re- 
ceiver output  through  a Federal  Scientific  hybrid  spectrum  analyzer 
and  outputting  the  results  on  an  ITT  facsimile  recorder.  The  analyzer 
applies  Hanning0®  weights  to  the  incoming  data  and  processes  a 0-to- 
2-kHz  bandwidth  to  4-Hz  resolution.  Figure  1H  shows  two  short 
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sequences  recorded  this  way.  Receiver  IF  frequency  is  plotted  versus 
time  (and  hence  azimuth)  with  amplitude  coded  as  intensity.  Vertical 
lines  are  indicative  of  strong  bursts  of  white  noise — presumably  from 
lightning.  Horizontal  lines  would  be  indicative  of  spectrally  pure 
signals  indicating  the  presence  of  some  types  of  interference. 

The  records  of  Figure  18  show  two  five-minute  periods  spaced  two 
hours  apart.  During  each  period  the  antenna  was  stepped  across  the 
same  8°  using  scanning  Format  A of  Table  1.  The  spectral  display  was 
useful  for  estimating  the  relative  amount  of  thunderstorm  activity 
versus  antenna  bearing.  For  example,  in  Figure  18(a),  the  northern 
bearings  are  the  most  active.  In  contrast,  two  hours  later 
[Figure  18(b)]  the  southern  bearings  have  become  extremely  active. 

As  we  shall  see,  this  later  activity  correlates  in  time  and  direction 
to  a known  severe  thunderstorm.  But  as  previously  described,  the 
observation  of  high  atmospheric  noise  activity  alone  cannot  provide  a 
reliable  indication  of  remote  severe  thunderstorms.  For  this  reason 
the  collected  data  were  next  analyzed  for  a class  of  HF  atmospherics 
that  might  be  uniquely  characteristic  of  severe  thunderstorms. 

C . Initial  Observations  of  an  Identifiable  Class  of  HF  Atmospherics 

Associated  with  Severe  Thunderstorms 

Each  day  the  National  Severe  Storms  Forecast  Center  compiles  a log 
noting  the  positions  and  times  of  all  severe  storm  sightings.  These 
logs  were  obtained  for  days  when  experimental  data  had  been  collected. 
Initial  analysis  consisted  of  reviewing  data  recorded  in  the  direction 
of  sighted  tornadoes.  The  ampli tude-versus-time  waveform  was  carefully 
reviewed  on  an  oscilloscope.  As  with  the  real-time  spectral  records, 
this  review  could  easily  discern  the  relative  level  of  storm  activity 
versus  antenna  bearing.  More  importantly,  at  times  and  bearings  coin- 
cident with  reported  severe  storm  manifestations  an  unusual  amount  of 
isolated  noise  impulses  was  observed.  This  was  the  first  indication 
of  a class  of  atmospherics  having  a remotely  observable  characteristic 
that  could  be  uniquely  correlated  with  tornadoes  and  possibly  other 
remote  severe  storm  events. 


As  described  in  Section  I,  HF  noise  from  lightning  is  character- 
ized by  many  impulses  (each  lasting  only  microseconds)  radiating  for 
periods  of  ten  to  hundreds  of  milliseconds.  For  communication  receiver 
bandwidths  (such  as  those  used  in  this  experiment),  individual  impulses 
are  not  resolved  and  the  atmospheric  appears  as  a quasi-continuous 
"burst"  of  noise.  This  is  the  type  of  in-beam  atmospheric  commonly 
observed  with  the  WARF  antenna.  While  the  time/ amplitude  details  of 
each  burst  differ,  they  seldom  last  less  than  10  ms.  Thus  the  observa- 
tion of  many  shorter  impulses  cannot  be  explained  on  the  basis  of 
normal,  even  active  thunderstorm  activity. 

Initially,  data  collected  on  19  May  1975  were  digitized  for  fur- 
ther analysis.  Several  10-s  records  from  these  data  are  plotted  in 
Figure  19.  In  each  record,  1-ms  samples  of  received  noise  magnitude 
are  plotted  versus  time.  One  second  of  data  appears  on  each  line  with 
subsequent  seconds  plotted  one  above  the  next.  The  three  10-s  records 
of  Figure  19  were  taken  at  the  times  and  bearings  noted.  Figure  19(b) 
was  recorded  at  the  time  and  in  the  direction  of  a tornado  sighted  near 
San  Antonio,  Texas  (a  distance  of  2200  km  from  the  WARF  antenna).  The 
other  two  records  were  recorded  at  times  or  in  directions  not  associ- 
ated with  reported  tornado  activity.  In  all  three  records,  atmo- 
spherics characteristic  of  normal  lightning  can  be  observed.  Superim- 
posed on  this  background  the  record  of  Figure  19(b)  contains  a large 
number  of  isolated  impulses.  While  Figures  19(a)  and  (c)  contain  some 
impulses,  the  number  and  amplitudes  of  those  in  19(b)  arc  significantly 
higher.  Based  on  this  observation,  a study  was  conducted  to  determine 
the  impulse  content  of  the  recorded  noise, 

I) , A Computer  Algorithm  for  Counting  Isolated  Impulses 

To  measure  the  impulse  content  of  atmospheric  noise  on  a large 
data  set  it  was  necessary  to  devise  some  automated  impulse-counting 
algorithm.  An  algorithm  employing  a fixed  threshold  would  count  all 
noise  bursts  and  the  measurement  would  be  dominated  by  the  longer  atmo- 
spherics possibly  originating  from  normal  in-beam  thunderstorms.  Thus, 
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any  counting  algorithm  developed  must  be  able  to  count  isolated  im- 
pulses (presumed  to  originate  uniquely  within  severe  storms)  while 
ignoring  those  bursts  lasting  several  milliseconds  or  more  (likely  to 
originate  within  normal  storms).  By  normalization  of  each  1-ms  data 
sample  by  a short-term  (10  to  50  ms)  estimate  of  its  surrounding  noise 
level,  a constant  threshold  algorithm  was  devised  to  count  only  iso- 
lated impulses. 

Strong  noise  samples  from  normal  lightning  flashes  are  surrounded 
in  time  by  other  strong  noise  samples.  For  this  reason  the  short-term 
average  noise  level  is  high  and  a sample  normalized  by  this  average  is 
not  likely  to  exceed  a preselected  threshold.  For  isolated  impulses, 
the  surrounding  average  noise  level  is  low  and  the  normalized  sample 
is  likely  to  be  counted. 

Figure  20  shows  the  effect  of  this  impulse-counting  algorithm.  A 
1— s record  of  noise  is  shown  in  Figure  20(a).  Figure  20(b)  plots  the 
short-term  (20  ms)  running  noise  average.  Figure  20(c)  replots  the 
noise  of  Figure  20(a)  but  normalized  by  the  appropriate  value  of  20(b). 
Also  drawn  is  a constant  12-dB  threshold.  Note  that  the  only  points 
to  exceed  the  threshold  are  from  isolated  impulses.  The  algorithm 
effectively  discriminated  against  the  strong  and  prolonged  atmospheric 
at  the  end  of  this  record. 

The  data  of  19  May  1975  were  used  to  test  various  choices  of 
averaging  time  and  threshold  in  an  attempt  to  optimize  the  algorithm's 
immunity  to  bursts  while  remaining  sensitive  to  impulses.  Averaging 
times  between  10  and  50  ms  worked  well,  while  shorter  times  contained 
too  many  spurious  counts  from  normal  lightning,  and  longer  times  tended 
to  miss  too  many  isolated  impulses.  Threshold  values  between  0 and 
20  dB  were  also  tried.  As  the  threshold  was  increased,  the  number  ol 
impulses  counted  decreased.  However,  all  such  counts  decreased  pro- 
portionally, leaving  the  relative  count  of  impulses  roughly  unchanged. 
Thus  the  algorithm  proved  insensitive  to  choices  of  parameters  over  a 
substantial  range  of  values.  A 20-ms  averaging  period  and  a 12-dB 
threshold  were  selected  for  analyzing  the  full  data  set.  Results  are 
presented  in  the  next  section. 
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FIGURE  20  FEATURES  OF  THE  IMPULSE -COUNTING  ALGORITHM 


VI  HF  IMPULSIVE  NOISE  CORRELATED  WITH 
REMOTE  SEVERE  THUNDERSTORMS 


In  this  section,  results  of  the  impulse-counting  algorithm  are 
compared  with  reported  sightings  of  severe  thunderstorm  events.  Data 
collected  on  the  first  day  of  the  experiment,  19  May  1975,  are  pre- 
sented in  detail.  Test  results  from  other  days  are  also  summarized. 
The  outcome  of  a statistical  test  performed  to  verify  the  observed 
correlation  is  described.  Finally,  the  validity  of  procedures  and 
assumptions  used  to  formulate  the  statistical  test  is  discussed. 

A.  Analysis  of  Data  Collected  on  19  May  1975 

The  numbered  dots  in  Figure  21  indicate  locations  of  tornadoes 
sighted  in  the  vicinity  of  San  Antonio,  Texas  on  19  May  1975.  Data 
collected  at  WARF  on  that  afternoon  were  processed  by  the  counting 
algorithm  described  in  the  previous  section.  The  results  are  plotted 
in  Figure  22.  Pulse  counts  are  plotted  versus  time  for  each  of  the 
five  antenna  bearings.  Times  of  the  sighted  tornadoes  are  indicated 
on  the  baseline  of  the  closest  bearing.  The  numbers  correspond  to 
those  of  Figure  21.  A strong  correlation  between  high  impulse  counts 
and  tornado  sightings  is  evident.  In  general,  as  the  storm  activity 
increased,  the  impulse  count  increased,  reaching  local  maxima  near  the 
sighting  times.  At  the  time  these  data  were  collected,  this  unusual 
impulsive  noise  component  had  not  yet  been  associated  with  tornado 
activity.  In  the  absence  of  current  weather  information,  data  collec- 
tion was  unfortunately  terminated  during  the  height  of  storm  activity. 

During  the  period  from  approximately  2225  UT  until  2300  UT,  high 
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impulse  counts  were  observed  from  the  bearing  of  102  T.  No  tornado 
activity  was  sighted  at  that  bearing.  Three  possible  explanations  can 
be  proposed  for  this  conflicting  observation: 

(I)  These  nigh-impulse  counts  represent  a ''false  alarm"  caused 
by  random  noise,  interference,  or  atmospherics  from  normal 
storms . 
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FIGURE  21  WARF  ANTENNA  SCAN  SECTOR  IN  RELATION  TO  TEXAS  SEVERE  STORM 
SIGHTINGS 


The  WARF  observations  correctly  indicated  the  presence  of 
a severe  storm,  but  no  observable  manifestations  were 
sighted  in  the  vicinity  of  the  storm. 

The  impulses  originated  from  the  severe  thunderstorm  that 
spawned  the  tornado  at  Location  1 in  Figures  21  and  22 
(matched  in  time  but  off  4°  in  azimuth)  and  were  received 
via  a close-i  i side lobe. 
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FIGURE  22  A COMPARISON  OF  IMPULSE  COUNTS  WITH  SEVERE  STORM 
SIGHTINGS  — 19  MAY  1975 
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Analysis  oi'  the  temporal  and  spectral  details  of  these  data  shows 
no  evidence  of  interference.  While  it  is  generally  held  that  tornadoes 
often  pass  unnoticed,  it  would  be  counterproductive  to  excuse  all  data 
inconsistencies  by  Explanation  2,  above.  Results  from  other  data  plus 
the  subsequent  discovery  of  a minor  antenna  cabling  fault  strongly 
suggest  that  atmospherics  from  the  storm  that  produced  the  tornado  at 
Location  1 were  received  via  a close-in  antenna  sidelobe. 

The  initial  qualitative  review  of  the  data  from  19  May  1975  estab- 
lished a significant  correlation  between  an  identifiable  class  of  atmo- 
spherics and  known  remote  tornadoes.  In  order  to  substantiate  this, 
correlation  data  recorded  on  three  other  days  were  similarly  analyzed. 


B. 


Results  of  the  Counting  Process 


A total  of  23  hours  of  data  collected  on  four  days  (including 

19  May  1975)  were  analyzed.  The  resulting  impulse  counts  are  plotted 

in  Figures  22  through  25.  Sightings  of  severe  weather  reported  either 
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by  the  National  Severe  Storm  Center  or  Storm  Data  are  indicated  by 
T for  tornado,  F for  funnel,  H for  hail,  and  G for  wind  gusts.  Specific 
comments  concerning  these  figures  are  given  in  the  following: 

(1)  Figure  23 — 22  May  1975.  Although  not  in  regions  scanned  by 
the  WARF  antenna,  many  tornadoes  occurred  during  the  record- 
ing period.  The  lack  of  observations  of  high  impulse  counts 
suggests  that  the  antenna  sidelobe  levels  remained  good  for 
bearings  well  removed  from  the  main  lobe.  The  high  impulse 
count  at  2040  UT  and  74  T was  caused  by  a momentary  burst  of 
impulsive  interference.  In  general,  the  counting  algorithm 
was  insensitive  to  low-level  interference  of  most  types. 

(2)  Figure  24 — 23  May  1975 . While  good  temporal  correlation  is 
evident,  again  there  are  indications  of  close-in  antenna 
sidelobes . 

(3)  Figure  25 — 28  May  1975.  On  this  day,  data  were  collected 
for  30  s at  each  of  ten  bearings  (Format  B,  Table  1). 

While  the  correlations  observed  in  these  plots  are  not  as  striking  as 
those  from  19  May,  in  general,  high  impulse  counts  correlate  in  time 
and  approximately  in  azimuth  with  severe  storm  sightings. 


14  4 MHz 


8 


O 


FIGURE  23  A COMPARISON  OF  IMPULSE  COUNTS  WITH  SEVERE  STORM  SIGHTINGS  — 22  MAY  1975 
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C , Testing  the  Statistical  Significance  of  the  Observed  Correlations 

The  plots  of  Figures  22  through  25  show  a general  correlation  be- 
tween high  impulse  counts  and  remote  severe  thunderstorm  events.  A 
more  reliable  assessment  of  this  correlation  can  be  obtained  through 
statistical  analysis. 

The  term  correlation  has  been  loosely  used  to  describe  a relation 
between  reception  of  high  counts  of  impulses  and  the  occurrence  of  mid- 
western  severe  thunderstorms.  If  we  consider  the  observed  impulse 
counts  to  be  caused  by  a random  process,  we  can  precisely  define  cor- 
relation as  a statistically  meaningful  difference  between  the  proba- 
bility density  function  governing  the  impulse  counts  observed  during 
periods  of  normal  storm  activity  and  that  governing  the  impulse  counts 
observed  during  periods  of  severe  storm  activity.  Specifically,  im- 
pulse counts  observed  during  severe  thunderstorms  have  a higher  mean. 

To  prove  this  we  must  test  (and  reject)  the  hypothesis  that  the  two 
random  processes  are  identical. 

1 . Parti tioning 

To  test  for  the  equality  (or  difference)  between  the  random 
processes  generating  impulse  counts  during  periods  with  or  without 
severe  storm  activity,  the  observations  must  first  be  partitioned  into 
these  two  categories.  Sightings  are  the  only  source  of  independent 
information  available  for  defining  these  partitions.  Unfortunately, 
severe  storm  events  often  go  unnoticed.  Even  when  observed,  sightings 
provide  incomplete  information  on  the  duration  and  location  of  each 
event.  For  these  reasons  partitioning  errors  cannot  be  avoided.  Par- 
titioning procedures  were  selected  to  maintain  a sensitive  indication 
of  storm  activity,  while  introducing  as  little  bias  as  possible. 

First,  to  simplify  data  analysis,  all  typos  of  severe  events 
(o.g.,  tornadoes,  hail,  etc.)  were  lumped  together.  Second,  azimuth 
information  was  discarded  to  avoid  partitioning  difficulties  caused  by 
inexact  knowledge  of  the  spatial  position  of  the  storm  and  the  apparent 
close-in  sidelobcs.  Azimuth  information  was  discarded  by  selecting  the 
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maximum  impulse  count  within  each  azimuth  scan  to  represent  that  scan. 
Collapsing  the  data  in  this  manner  has  the  effect  of  making  time  of 
day  the  only  correlation  parameter. 

To  partition  the  data,  a time  interval  corresponding  to  each 
sighting  must  be  selected.  Because  the  duration  of  events  can  differ, 
selection  of  a fixed  partitioning  interval  is  inexact,  but  the  lack  of 
detailed  supporting  information  leaves  no  alternative.  A partitioning 
interval  of  ±10  minutes  centered  about  the  sighting  times  was  selected 

With  the  use  of  these  partitioning  rules,  the  data  from 
Figures  22  through  25  were  partitioned  into  ''normal''  and  "severe"  cate- 
gories. Histograms  of  the  resulting  observations  were  plotted  in 
Figure  26.  These  histograms  show  marked  differences.  The  statistical 
procedure  described  below  was  used  to  test  against  the  possibility 
that  these  differences  could  have  happened  merely  by  chance. 


FROM  NORMAL  PARTITION  (160  observations) 
FROM  SEVERE  PARTITION  (124  observations) 
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FIGURE  26  HISTOGRAMS  OF  IMPULSE  COUNTS  FOR  SEVERE  AND  NORMAL  PARTITIONS 
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2. 


Statistical  Hypothesis  Testing 


Let  the  impulse  counts  partitioned  into  the  normal  category 
be  labeled  X^  ...  X . They  shall  be  considered  samples  from  a popula- 
tion, X,  with  probability  density  function  f^fZ).  Similarly,  let  the 
impulse  counts  partitioned  into  the  severe  category  be  labeled  ...  Y. 
They  shall  be  considered  samples  from  a population,  Y,  with  probability 

density  function  f (Z).  The  null  hypothesis,  H , can  then  be  stated  as 
Y o 

follows.  The  probability  density  functions  of  the  normal  and  severe 
populations  are  identical — i.e., 

fx(Z)  = fy(Z)  for  all  Z.  (1) 

The  validity  of  the  null  hypothesis  can  be  tested  against  the  alterna- 
tive, H^,  which  is  defined  as  follows.  The  probability  density  func- 
tions of  the  severe  population  is  greater  than  that  of  the  normal 
population  (for  at  least  one  value  of  Z) — i.e., 

f (Z)  < f (Z)  for  at  least  one  Z.  (2) 

X.  Y 


The  null  hypothesis  is  accepted  or  rejected  by  comparing  a 
test  statistic,  U,  (based  on  the  observations)  to  a threshold,  T (based 
on  the  appropriate  probability  distributions  and  the  desired  level  of 
significance).  The  level  of  significance,  a,  determines  the  probabil- 
ity of  falsely  rejecting  11^ . For  the  tests  described  below,  an  a of 
0.01  was  chosen.  This  implies  that  if  there  were  truly  no  correlation, 
the  chances  are  one  in  one  hundred  of  erroneously  concluding  that  there 
was  a correlation. 

If  the  probability  density  functions  (1\.  and  f ) are  known, 
a parametric  test  can  be  used.  Since  we  wish  to  minimize  the  number  of 
assumptions  made,  a non-parametric  test  or  the  hypothesis  was  chosen. 
While  such  tests  are  slightly  less  exact  than  their  parametric  counter- 
parts, they  are  less  dependent  on  assumptions  and  usually  easier  to 
compute . 
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3 . A Non-Paranietric  Test  for  Correlation 

52 

The  Wilcoxon  Rank  sum  test  was  selected  as  an  appropriate 
non-parametric  test  of  the  stated  hypothesis.  To  compute  the  necessary 
test  statistic,  the  n observations  of  X and  the  m observations  of  Y are 
first  lumped  together  and  then  ordered.  A rank  is  assigned  to  each 
observation  according  to  its  position  in  the  order  (i.e.,  1 to  the 
smallest  observation  and  m + n to  the  largest).  The  sum  of  ranks,  T^., 
assigned  to  the  severe  category  is  next  computed.  Under  the  null  hypo- 
thesis the  two  probability  density  functions  are  equal  and  the  ranked 

observations  should  be  thoroughly  mixed.  Thus,  under  11  it  can  be 
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shown  “ that: 


E[Ty]  = 


m(m  + n -h  1) 


Var[TY]  = 


inn  (m  + n + 1 ) 


For  large  values  of  m and  n the  distribution  of  the  test  statistic, 


ty  - e[ty] 

/Var[  T "I 


is  closely  approximated  by  a standard  normal  distribution.  To  test 
H)(  the  computed  U is  compared  to  a threshold  value  determined  from 
the  0.01  significance  level  for  a standard  normal  distribution  function. 
Applying  this  test  to  the  partitioned  data  resulted  in  rejection  of  11  . 
Calculations  used  in  the  Wilcoxon  test  are  summarized  in  Table  2.  Re- 


jecting the  null  hypothesis  supports  the  conclusion  that  a statistically 
significant  correlation  exists  between  the  observed  impulse  counts  and 
the  remote  severe  storm  events.  It  should  be  emphasized  that  the 
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Tabic  2 

A SUMMARY  OF  COMPUTATIONS  FOR 
THE  W1LCOXON  RANK  SUM  TEST 


Number  of  counts  in  normal  category,  n 

160 

Number  of  counts  in  severe  category,  m 

124 

Sum  of  severe  category  ranks,  T^ 

24,133 

Expected  mean  of  Ty,  E[ T^l 

17,670 

Expected  standard  deviation  of  Ty, 

/Var[T  "! 

686.  M 

Test  statistic,  U 

9.  12 

Threshold  at  0.01  significance  level,  T 

2.33 

U = 9 .42  > T - 2.32 

Conclusion:  Reject  11 

() 

calculated  value  of  U is  very  high — more  than  nine  standard  deviations 
above  the  mean  of  a normal  distribution.  The  random  chance  of  calcu- 
lating such  a high  value  of  U when  H is  true  is  negligible. 


D . Defense  of  Assumptions 

In  this  section  results  were  presented  from  an  experiment  to 
determine  if  severe  storm  activity  radiates  an  identifiable  class  of 
HF  atmospherics  that  could  be  remotely  detected.  Statistical  tests 
were  used  to  show  that  the  observations  were  significantly  correlated 
with  known  severe  storm  events.  It:  is  necessary  to  review  these  re- 
sults in  light  of  the  assumptions  made  and  the  possibilities  of  bias. 

The  Wilcoxon  Rank  sum  test  assumes  that  the  two  testing  categories 
are  sampled  from  independent  populatin'  s . This  assumption  was  clearly 
violated  because  the  lack  of  exact  control  information  made  partition- 
ing errors  unavoidable.  However,  contamination  caused  by  partitioning 
errors  serves  only  to  mix  the  two  categories  and  hence  support  the  null 
hypothesi s . 


Besides  testing  the  validity  of  assumptions  we  must  be  careful  to 
ensure  that  no  systematic  bias  has  been  introduced  into  the  tests.  The 
method  of  partitioning  is  the  major  possible  source  of  bias.  The  col- 
lapsing of  azimuth  information  is  one  potential  source.  This  procedure 
was  selected  because  of  known  physical  uncertainties  in  the  data.  Fur- 
thermore it  can  be  argued  that,  under  the  null  hypothesis,  azimuth  col- 
lapsing should  not  effect  the  outcome.  A second  source  of  possible 
bias  comes  from  the  selection  of  a partitioning  interval.  But  as 
pointed  out  above,  the  use  of  a fixed  partitioning  interval  is  likely 
to  mix  the  partitions  and  hence  produce  a conservative  bias — i.e.,  one 
favoring  the  null  hypothesis. 

Based  on  the  foregoing  discussion  the  author  considers  the  statis- 
tical results  to  be  valid.  It  can  be  concluded  from  the  data  and  tests 
presented  here  that  the  impulsive  class  of  HF  atmospherics  is  radiated 
during  the  occurrence  of  severe  thunderstorm  events  such  as  tornadoes, 
and  that  this  noise  can  be  remotely  detected  on  a directional  antenna. 
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PROPERTIES  OF  HF  ATMOSPHERICS  RADIATED 
FROM  SEVERE  THUNDERSTORMS 

The  previous  two  sections  established  that  HF  atmospherics  of  a 
unique  class  were  radiated  during  severe  thunderstorms . In  this  sec- 
tion a comparison  between  these  atmospherics  and  results  from  similar 
line-of-sight  measurements  is  presented.  Results  of  this  comparison 
show  that  the  remote  observations  represent  a class  of  HF  atmospherics 
not  previously  identified.  The  temporal  characteristics  ol  this  new 
class  of  atmospherics  are  analyzed,  and  some  speculation  is  presented 
as  to  the  meteorological  processes  within  severe  thunderstorms  that 
might  generate  these  atmospherics. 

A . Comparison  of  Remote  and  Linc-of-Sight  HF  Atmospherics — 

Identification  of  a New  Class  of  Atmosphcri c s 

Line-of-sight  observations  of  HF  atmospherics  from  severe  thunder- 
storms have  been  made  by  Stanford, Lind,'^  and  Taylor  * In  their 

experiments  the  total  number  of  received  atmospheric  impulses  were 
counted.  Since  HF  atmospherics  from  normal  lightning  consist  of  hun- 
dreds of  impulses,  it  is  clear  why  counts  of  received  impulses  made 
by  these  researchers  were  typically  very  large.  Stanford  and  Lind  used 
this  count  directly  to  correlate  with  local  severe  storm  activity. 

Taylor  used  the  number  of  impulse  bursts  whose  rate  exceeded  a selected 
threshold . 

It  was  explained  in  Section  V why,  for  remote  measurements,  the 
total  number  of  received  atmospheric  impulses  is  not  a sensitive  mea- 
sure of  the  severity  of  any  single  storm.  Because  it  is  not,  the  ex- 
periment described  in  Section  V was  designed  to  search  for  a class  of 
severe  storm  atmospherics  that  could  be  extracted  from  the  bulk  of 
extraneous  in-beam  atmospherics.  The  class  of  atmospherics  subsequently 
discovered  and  correlated  with  severe  storms  can  be  described  as  iso- 
lated impulses  and  represent  only  about  271  of  the  total  atmospheric 
content  received. 
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Consequently,  the  isolated  impulse-counting  algorithm  was  designed 
to  reject  most  of  the  impulses  that  would  be  counted  by  Stanford  and 
Lind,  and  to  reject  all  of  the  impulse  bursts  measured  by  Taylor.  Con- 
versely, the  remotely  measured  isolated  impulses  would  account  for  only 
a very  small  fraction  of  the  impulses  measured  by  Stanford  and  Lind, 
and  would  never  exceed  the  burst-rate  criteria  of  Taylor.  On  the  basis 
of  these  considerations,  it  can  be  assumed  that  the  class  of  remotely 
received  HF  atmospherics  described  here  has  not  been  previously  identi- 
fied through  line-of-sight  measurements.  It  would  be  useful  to  apply 
these  different  measurement  procedures  to  the  same  line-ol'-sight  data 
to  determine  their  relative  sensitivity  to  severe  thunderstorms. 

B.  Temporal  Characteristics  of  the  Observed  Severe  Storm  Atmospherics 

Since  the  severe  storm  HF  atmospherics  of  the  type  described  here 
have  not  been  studied  before,  it  is  worthwhile  measuring  some  of  their 
characteristics.  While  it  would  be  desirable  to  study  these  atmo- 
spherics in  terms  of  their  frequency  spectrum,  amplitude  distribution, 
detailed  waveform,  etc.,  the  available  data  limit  analyses  made  thus 
far  to  the  time  sequence  of  the  impulses.  Still,  this  temporal  char- 
acterization can  be  useful  for  estimating  what  processes  might  cause 
such  unusual  atmospherics.  Accordingly,  the  time  sequence  ol  impulses 
was  analyzed  to  determine  if  it  followed  a Poisson  Law,  as  might  be 
expected  for  a completely  random  process,  or  if  it  showed  a more 
structured  distribution  as  might  be  expected  if  the  causal  process  was 

periodic,  as  some  line-of-sight  atmospherics  measurements  have  indi- 
27 

cated  in  the  past. 

It  follows  that  if  the  atmospherics  are  distributed  along  the 

time  axis  according  to  the  Poisson  Law,  then  the  time  intervals  between 
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atmospherics  arc  exponentially  distributed.  For  this  reason  the  time 
intervals  between  atmospherics  were  analyzed  by  plotting  their  cumula- 
tive distribution.  Figure  27  plots  this  distribution  for  a one-minute 
sample  recorded  during  the  19  May  1975  severe  thunderstorm  period. 

Also  plotted  is  an  idealized  exponential  distribution.  These  plots 


(5(t 


DATA 


2 3 

NORMALIZED  TIME  INTERVALS 


658  D 3?  961 


FIGURE  27  A CUMULATIVE  PROBABILITY  DISTRIBUTION  OF  THE  INTERVALS  BETWEEN 
IMPULSES  SUPERIMPOSED  ON  AN  EXPONENTIAL  DISTRIBUTION 


show  reasonable  agreement,  suggesting  that  the  atmospherics  are 
generated  by  a Poisson  process. 

But  before  these  atmospherics  can  be  modeled  as  a simple  Poisson 
process,  it  is  necessary  to  establish  that  they  are  independent  events 
and  that  their  statistics  do  not  vary  with  time  (i.c.,  are  stationary) 
This  latter  assumption  cannot,  of  course,  be  strictly  true  because,  as 
was  shown  in  the  last  section,  there  is  a definite  correlation  between 
the  number  of  atmospherics  received  and  the  occurrence  of  severe  thun- 
derstorm events.  However  these  storm  events  occur  on  a time  scale  ot 
tens  of  minutes,  and  it  would  be  useful  to  determine  whether  the  sta- 
tistics of  these  atmospherics  are  roughly  constant  for  short  periods 
(e.g.,  tens  of  seconds).  For  if  short-term  stationarity  can  be  estab- 
lished, the  atmospherics  can  lie  modeled  as  a Poisson  process  with  a 
slowly  varying  mean. 

To  test  for  independence  and  short-term  stationarity,  two  further 
calculations  were  performed  on  the  1ll  May  data.  First,  the  autocor- 
relation function  was  computed  from  the'  intervals  between  events  (see 


FIGURE  28  THE  AUTOCORRELATION  FUNCTION  OF  THE  TIME  INTERVALS  BETWEEN 
SEVERE  STORM  ATMOSPHERICS 


Figure  28).  For  an  independent  process  this  function  should  have  small 
values  for  all  lags  other  than  zero.  Second,  the  cumulative  count  of 
atmospherics  was  calculated  as  a function  of  time  as  shown  in  Figure  29. 
For  a stationary  process  this  plot  should  approximate  a straight  line. 

The  three  types  of  temporal  analyses  described  above  were  performed 
on  10  one-minute  samples  collected  during  the  storm  of  19  May  1975.  All 
data  samples  gave  essentially  the  same  positive  results  as  those  de- 
picted in  Figures  27,  28,  and  29.  While  this  is  not  an  exhaustive  test, 
it  does  provide  strong  indications  that  the  time  sequence  of  remotely 
observed  severe  thunderstorm  atmospherics  may  be  generated  by  a Poisson 
process  having  a slowly  varying  mean  related  to  the  severe  storm's 
intensity . 


C . Meteorological  Processes  Likely  to  Cause  the  Observed  Atmospherics 

At  this  point  it  may  be  beneficial  to  speculate  about  processes 
within  severe  thunderstorms  that  might  cause  the  observed  atmospherics. 
If  the  source  of  those  atmospherics  can  be  identified,  future  atmo- 
spheric measurements  might  be  used  to  study  physical  processes  within 
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FIGURE  29  THE  CUMULATIVE  COUNT  OF  SEVERE  STORM 
ATMOSPHERICS  AS  A FUNCTION  OF  TIME 


the  storm.  First,  let  us  consider  what  can  be  deduced  from  the  mea- 
surements. The  atmospherics  were  observed  at  HF  and  thus  were  likciy 
produced  from  discharges  of  approximately  20  meters.  To  the  1-ms 
resolution  achievable  with  the  equipment  used,  these  atmospherics 
appear  as  simple  impulses,  each  independent  of  the  rest.  These  char- 
acteristics are  very  different  from  HF  atmospherics  of  normal  lightning 
(comprised  of  many  impulses  following  in  close  succession).  From  these 
observations  we  can  assume  that  the  lightning  causing  these  atmospherics 
in  severe  storms  is  a class  of  frequently  occurring  but  rather  minor 
discharges.  Two  plausible  generating  mechanisms  will  be  considered 
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next  . 


Severe  thunderstorms  generally  contain  greater  turbulence  and 
stronger1  electrification  than  do  normal  thunderstorms.  Because  of  the 
added  turbulence,  let  us  assume  that  the  usual  dipole  distribution  of 
electric  charge  within  the  thundercloud  might  be  redistributed.  Winds 
could  carry  smaller  pockets  of  charge  randomly  about.  When  two  op- 
positely charged  pockets  encounter  one  another  a single  relatively 
small  discharge  might  occur.  This  process  would  happen  throughout  the 
cloud,  creating  the  random  atmospherics  observed.  Under  this  hypothesis 
the  rate  of  atmospherics  might  provide  information  on  the  amount  of 
turbulence  and  electrification  present  in  the  severe  storm. 

Under  a second  hypothesis  the  observed  atmospherics  would  be  caused 
by  a process  that  generates  its  own  free  electricity.  Recent  Doppler 
radar  measurements  have  identified  concentrated  regions  within  the 
storm  cloud  having  strong  wind  shears  and  vortex  motion.  These  local- 
ized winds  can  occur  10  to  20  minutes  before  a tornado  sighting.  in 
some  instances  tornadoes  have  been  observed  to  originate  from  these 
vortices.  It  is  possible  that  these  strong  winds  might  produce  their 
own  electricity.  If  they  do,  this  could  account  for  the  short  random 
atmospherics  observed.  The  rate  of  atmospherics  might  then  provide 
information  on  the  occurrence  and  intensity  of  these  localized  events. 

Of  the  two  hypotheses  presented  above,  1 he  author  considers  the 
latter  the  most  plausible.  This  Is  based  on  the  following  considera- 
tions. First,  the  turbulence  theory  would  suggest  that  those  atmo- 
spherics might  he  routinely  observed  in  all  storms.  While  I lie  data 
presented  here  arc  inconclusive,  it  is  the  author’s  opinion  that  these 
atmospherics  are  not  consistently  present  during  normal  storms.  Second, 
the  turbulence  theory  assumes  a general  process  occurring  throughout  tin 
storm  cloud.  Thus,  the  atmospherics  would  lit'  proportional  to  tin  over- 
all storm  intensity.  However,  the  data  show  strong  enhancements  in 
the  number  of  received  atmospherics  occurring  at  the  time  of  sighted 
severe  storm  events  (e.g.,  see  Figure  22).  The  observ  >tions  therefore 
are  more  consistent  with  the  winel-shear  theory,  which  suggests  that 
these  atmospherics  would  be  more  closely  linked  with  specilic  and 
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localized  storm  events.  Lastly,  there  have  been  observations  of 

24 

lightning  discharges  within  the  actual  tornado  funnel.  Since  the 

funnel  can  be  considered  as  an  extension  of  the  in-cloud  vortex,  its 

lightning  might  also  be  caused  by  extreme  wind  shears. 

Section  IV  discussed  an  HF  direction-finding  technique  that 

demonstrated  that  the  distribution  of  lightning  within  a cloud  could 

be  remotely  measured.  After  suitable  modification  for  processing  the 

shorter  severe  storm  atmospherics,  this  technique  could  be  used  to 

measure  the  angular  distribution  of  these  atmospherics.  If  they  were 

found  to  be  distributed  over  a wide  region,  this  would  support  the 

hypothesis  that  they  are  caused  by  a general  increase  in  turbulence. 

On  the  other  hand,  if  they  were  found  to  be  tightly  clustered,  this 

* 

would  support  the  wind-shear  hypothesis. 

In  concluding  this  section  it  must  be  emphasized  that  the  above 
discussion  is  largely  speculative.  It  is  included  to  present  the 
reader  with  some  unproven  thoughts  and  to  demonstrate  that  the  remote 
observation  of  these  atmospherics  can  be  useful  for  studying  physical 
processes  occurring  within  the  storm  cloud. 
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It  should  be  noted 
have  been  conducted 
not  likely  that  VLK 
crates  the  observed 


that  line-of-sight  direction-finding  measurements 
at  VLF,  with  inconclusive  results.  However,  it  i 
radiation  is  caused  by  the  same  process  that  gen- 
HF  atmospherics. 
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VIII  POSSIBILITIES  FOR  THE  REMOTE  STUDY  OF 
LIGHTNING  AND  RELATED  METEOROLOGICAL  PROCESSES 


The  research  results  described  in  this  dissertation  provide  new 
and  hopefully  useful  information  about  severe-thunderstorm  electricity. 
This  research  was  conducted  because  it  was  recognized  that  the  proper- 
ties of  a directional  HF  antenna  were  well  suited  for  remote  measure- 
ments of  severe  storm  atmospherics.  In  this  section  we  shall  briefly 
consider  how  the  capabilities  of  a directional  HF  antenna  can  be  ap- 
plied to  gather  further  information  that  is  not  conveniently  obtainable 
by  other  measurement  techniques. 

The  most  useful  remote  sensing  capability  of  a directional  antenna 

is  the  ability  to  select  which  storms  within  a very  wide  area  will  be 

monitored.  For  example,  the  WARF  antenna  can  receive  one-hop  atmo- 
* 

spherics  from  storms  within  a region  from  Texas  to  South  Dakota  and 
the  Rocky  Mountains  to  the  Mississippi  River.  Since  most  U.S.  torna- 
does occur  within  this  region,  the  remote-sensing  technique  can  provide 
a relatively  inexpensive  method  to  observe  a large  number  of  geograph- 
ically separated  severe  thunderstorms  within  a single  storm  season. 

Such  a measurement  program  could  determine  the  percentage  of  severe 
storms  exhibiting  identifiable  electrical  activity  and  hence  determine 
the  value  of  HF  atmospherics  for  tornado  warnings. 

The  ability  to  monitor  a large  region  is  a second  attractive  capa- 
bility of  this  remote-measurement  technique,  because  wide-area  measure- 
ments can  help  in  finding  and  continuously  monitoring  severe  thunder- 
storms. Wide-area  measurements  might  also  allow  for  the  study  of  the 
occurrence  rates  of  severe  thunderstorms  versus  geographical  location 
or  meteorological  conditions.  Finally,  in  contrast  with  line-of-sight 
measurements,  for  the  remote  measurements  the  movements  of  a storm 
cause  only  a small  percentage  change  in  range  and  aspect  angle.  These 
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This  region  is  substantially  reduced  during  times  when  blanketing 
sporadic-E  propagation  is  present. 
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factors  can  allow  for  the  observation  of  a thunderstorm  during  the 
entire  period  it  is  radiating  identifiable  atmospherics,  independent 
of  the  distance  the  storm  may  travel. 

For  any  single  storm,  normal  summertime  ionospheric  propagation 

* 

will  limit  the  usable  frequency  band  to  a few  megahertz.  However, 
this  may  be  sufficient  to  gain  a crude  estimate  of  the  lightning  dis- 
charge length.  This,  in  turn,  could  pi-ovidc  information  about  the 
average  discharge  length  versus,  for  example,  the  occurrence  of 
funnels . 

Lastly,  the  direction-finding  technique  described  in  Section  IV 
may  be  valuable  for  obtaining  a spatial  distribution  of  the  severe 
storm  atmospherics.  No  such  determination  of  the  location  of  dis- 
charges in  the  20-meter-wavelengtli  region  has  been  previously  at- 
tempted. It  was  pointed  out  in  Section  VII  how  such  measurements 
might  help  explain  the  mechanisms  causing  these  isolated  atmospherics. 

Thus  far,  I have  discussed  the  remote-sensing  properties  of  a 
directional  antenna  through  examples  having  applications  in  severe- 
thunderstorm  research.  There  may  also  be  several  uses  of  remote  mea- 
surements of  normal  storms.  For  example,  while  ample  line-of-sight 
data  have  been  collected  on  over-land  thunderstorms , virtually  no 
direct  measurements  have  been  made  on  open-ocean  thunderstorms.  And 
those  remote  observations  made  by  receiving  VLF  atmospherics  are 
largely  limited  to  the  return  stroke  of  a cloud- to-ground  flash.  How- 
ever, remote  measurements  of  HF  atmospherics  can  be  used  to  observe 
other  lightning  processes.  For  example,  atmospherics  from  both  cloud- 
to-ground  and  intra-cloud  flashes  can  be  received.  The  cloud- to-ground 
flashes  can  be  identified  by  observation  of  the  return  stroke  quenching 
period  described  in  Section  [V.  Thus,  the  proportion  of  over — ocean 
flashes  going  to  ground  can  bo  obtained.  Over-land  thunderstorms  show 
a relation  between  the  proportion  of  flashes  going  to  ground  and  the 
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Sporadic-E  propagation  or  winter-day  F- layer  propagation  could  provide 
as  much  as  a 2.5-to-l  usable  frequency  band. 
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storm's  latitude.  For  storms  in  the  WARF  coverage  area  (between 
25°N  and  45°N  altitude)  the  proportion  of  flashes  going  to  ground 
should  vary  from  17%  to  34%.  Testing  this  relationship  for  over-ocean 
storms  might  provide  insight  into  the  mechanisms  causing  this  lightning 
variation  with  latitude. 

In  summary,  it  must  be  emphasized  that  the  ability  to  accomplish 
the  measurements  proposed  above  has  not  been  fully  demonstrated.  How- 
ever, on  the  basis  of  the  measurements  that  have  been  accomplished,  the 
author  considers  these  further  applications  as  realistic  possibilities. 
It  is  suggested  that  the  remote  measurement  of  HF  atmospherics  by  a 
directional  antenna  can  be  a low-cost  method  of  studying  thunderstorm 
activity  for  a wide  range  of  scientific  and  applied  purposes. 
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IX  SUMMARY  AND  CONCLUSIONS 
A.  Summary  of  Results 

With  no  local  thunderstorms  present  it  was  observed  that  strong 
HF  atmospherics  from  individual  lightning  flashes  were  commonly  re- 
ceived on  a directional  HF  antenna.  In  contrast,  these  same  atmo- 
spherics were  only  very  weakly  received  on  a nearly  collocated  omni- 
directional antenna.  Furthermore,  the  amplitude-probability  distribu- 
tions calculated  from  noise  simultaneously  received  on  these  two 
antennas  showed  a significantly  stronger  high- amplitude  (and  impulsive) 
component  for  the  directional  antenna.  These  observations  suggested 
that  the  directivity  and  sidelobe  level  of  this  directional  antenna 
enabled  the  reception  of  atmospherics  from  one-hop  distant  in-beam 
lightning  flashes  with  sufficient  strength  above  background  noise  to 
permit  measurement  of  their  details. 

Direct  verification  of  this  hypothesis  was  obtained  during  an  ex- 
periment in  which  simultaneous  measurements  were  recorded  with  the 
directional  antenna  and  at  a field  site  situated  near  a remote  thunder- 
storm. The  directional  antenna  was  pointed  toward  the  field  site  (and 
hence,  the  thunderstorm),  which  was  at  a one-hop  distance.  The  simul- 
taneous recordings  showed  that  some  of  the  strong  atmospherics  received 
cn  the  directional  antenna  originated  from  the  remote  storm,  thus  mark- 
ing the  first  time  that  the  source  of  skywave  HI'  atmospherics  was  pre- 
cisely located.  Analysis  of  these  records  demonstrated  that,  despite 
undergoing  ionospheric  propagation,  these  skywave  atmospherics  were 
received  with  enough  strength  and  detail  to  provide  usable  information 
on  individual  remote  lightning  flashes. 

By  the  application  of  digital  beamforming  techniques  to  the  direc- 
tional antenna,  the  direction  of  arrival  of  individual  atmospherics  was 
measured  at  4-ms  intervals  with  0.06°  azimuthal  accuracy.  This  unusual 
measurement  revealed  that  the  direction  of  arrival  of  individual  atmo- 
spherics appeared  to  shift,  indicating  changes  in  the  horizontal  posi- 
tions of  in-cloud  radiators.  For  one  storm,  the  average  horizontal 
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dimension  of  the  electrically  active  region  involved  in  each  flash  was 
measured  to  be  9 ± 1 km.  This  is  comparable  to  the  3-to-15-km  range 
of 'values  previously  reported  by  line-of-sight  observations,  While 
there  is  no  direct  verification  that  the  observed  variations  in  direc- 
tion of  arrival  were  caused  by  changes  in  the  radiating  regions  within 
a single  flash,  calibration  signals  received  from  a transmitter  close 
to  the  thunderstorm  showed  that  these  unusual  observations  were  not 
equipmental  or  ionospherically  induced. 

The  directional  antenna  was  next  used  to  study  atmospherics  from 
remote  severe  thunderstorms.  An  experiment  was  conducted  to  determine 
if  midwestern  severe  thunderstorms  radiate  identifiable  atmospherics 
that  could  be  remotely  received.  Such  a class  of  atmospherics  was 
observed  during  periods  of  severe  storm  activity.  An  algorithm,  de- 
vised to  extract  this  class  of  atmospherics  from  other  noise,  was  ap- 
plied to  23  hours  of  data  collected  on  four  days  and  results  were  com- 
pared with  reports  of  severe  storm  sightings.  By  the  use  of  a non- 
parametric  statistical  test,  it  was  found  that  the  strong  presence  of 
this  class  of  atmospherics  correlated  remarkably  well  with  the  pres- 
ence of  severe  thunderstorm  events. 

These  observations  were  compared  with  published  results  of  HF 
measurements  made  at  linc-of-sight  distances  from  severe  storms.  It 
was  found  that  the  line-of-sight  processing  techniques  used  to  estimate 
severe  storm  activity  are  quite  different  from  the  technique  developed 
here.  On  the  basis  of  this  comparison,  it  is  not  likely  that  the  class 
of  remotely  observed  atmospherics  has  been  previously  isolated.  A 
study  of  the  temporal  characteristics  of  these  atmospherics  revealed 
that  they  closely  lit  a Poisson  model.  Finally,  it  was  speculated 
that  extreme  wind  shears  and  in-cloud  vortices  (and  perhaps  even  actual 
tornado  vortices)  found  in  severe  but  not  normal  thunderstorms  are  the 
likely  generating  mechanism  for  these  unusual  atmospherics . 
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Recommendations  for  Future  Work 


The  current  work  has  established  the  utility  of  a directional  HF 
antenna  for  remote  measurements  of  atmospherics  and  has  discovered  a 
new  class  of  atmospherics  associated  with  severe  thunderstorms.  Pos- 
sible further  applications  of  this  measurement  technique  have  already 
been  discussed.  Presented  below  are  specific  extensions  of  the  present 
work  applied  to  several  different  problem  areas. 

1 . The  Study  of  Severe  Thunderstorm  Atmospherics 

Further  remote  measurements  of  severe  storm  atmospherics 
should  be  directed  toward  obtaining  information  on  their  amplitude 
distribution,  frequency  spectrum,  location  within  the  storm,  and  the 
relation  of  these  characteristics  with  observed  severe  storm  events. 
Secondly,  an  attempt  should  be  made  to  measure  this  class  of  severe 
storm  atmospherics  at  line-of-sight  distances.  A comparison  of  the 
sensitivity  and  false-alarm  rates  of  the  different  algorithms  used  to 
evaluate  severe  storm  activity  would  be  useful. 

2 . The  Study  of  Atmospherics  from  Hurricane  Thunderstorms 

Little  is  kr  wn  about  the  thunderstorm  activity  involved 
with  off-shore  hurricanes.  Since  Gulf  of  Mexico  hurricanes  can  come 
within  one-hop  distance  of  the  WARF  antenna,  they  should  be  monitored 
for  any  indication  of  the  presence  of  unusual  HF  atmospherics. 

3 . The  Effects  of  Atmospherics  on  Adaptive  Boamfonning 

A topic  currently  receiving  much  attention  is  the  use  of 
adaptive  beamforming  to  increase  system  performance  on  large  HF  an- 
tennas. Further  characterization  of  the  ;unplitude,  temporal,  and  azi- 
muthal properties  of  atmospherics  received  on  large  antennas  would  be 
of  profit  to  this  technology. 
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Appendix  A 

THE  USE  OF  WIDE- SWEEP  BACKSCATTER  SOUNDINGS 
FOR  FREQUENCY  SELECTION 

During  the  past  decade  wide-sweep  backscatter  soundings  have  de- 
veloped into  a useful  ionospheric  diagnostic  technique.  Taken  from 
* 

one  location,  these  soundings  provide  information  on  current  propaga- 
tion conditions  for  frequencies  in  the  HF  band  and  for  distances  of 
interest.  A block  diagram  of  WARF  equipment  used  to  generate  these 
soundings  is  shown  in  Figure  A-l.  An  example  of  a wide-sweep  back- 
scatter sounding  is  reproduced  in  Figure  A-2 . For  these  soundings, 
backscattered  energy  is  coded  as  intensity  and  plotted  versus  radio 
frequency  and  propagation  time  delay.  Such  backscatter  soundings  were 
! periodically  taken  during  the  experiments  described  in  this  report. 

i 

Details  of  the  transmitted  waveform  and  data  processing  necessary  to 

5 6 

produce  such  soundings  have  been  described  by  Sweeney,  Barnum,  and 

I 7 

Washburn.  Of  concern  here  is  the  manner  in  which  these  soundings  are 
used  to  determine  an  appropriate  receiving  frequency. 

The  soundings  were  used  to  choose  a frequency  interval,  about 
1 MHz  wide,  that  would  provide  good  propagation  from  a desired  location 
or  region.  The  relative  propagation  strength  can  be  measured  from  the 
sounding  if  we  assume  the  scattering  cross  section  of  the  ground  to  be 
reasonably  constant.  This  is  an  excellent  assumption  for  sea  back- 
scatter and  remains  usable  for  land  backscatter  if  localized  enhance- 
ments caused  by  land  prominences  are  ignored. 

In  the  evaluation  of  propagation  conditions  for  a specific  loca- 
tion, the  ionospheric  reflecting  height  must  be  approximately  known. 
This  is  necessary  so  that  the  ground  range  to  the  location  can  be  con- 
verted to  propagation  time-delay  for  comparison  with  the  backscatter 
sounding.  An  estimate  of  the  ionospheric  reflection  height  can  be 


* 

Transmitter  and  receiver  may  be  slightly  separated  (o.g.,  separation 
for  the  WARF  sounder  is  approximately  185  km). 
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FIGURE  A-2  A WIDE-SWEEP  BACKSCATTER  SOUNDING 


obtained  by  identifying  the  propagation  layer  through  analysis  of  the 
wide-sweep  soundings  or  by  extrapolation  from  heights  measured  directly 
overhead . 


Once  the  height  of  reflection  has  been  estimated,  the  propagation 
time-delay  interval  for  a specific  ground-range  interval  (i.e.,  a 
specified  geographic  region)  can  be  determined.  The  sounding  is  used 
to  identify  the  frequency  interval  with  the  strongest  backscatter  over 
the  specified  time-delay  interval.  A radio  receiver  is  next  scanned 
across  the  selected  frequency  interval  to  choose  a single  channel  (of 
bandwidth  1.4  kHz)  that  is  free  of  interference. 

For  example,  the  sounding  of  Figure  A-2  was  taken  on  23  May  when 
the  weather  watch  of  Figure  17  was  in  effect.  Analysis  of  the  sounding 
suggested  that  the  principal  mode  of  ionospheric  propagation  was 
F-layer.  The  overlay  of  Figure  17  was  constructed  assuming  typical 
F-layer  heights  and  was  used  to  convert  the  geographic  location  of  the 
weather  watch  into  a time-delay  interval.  From  the  sounding  of 
Figure  A-2  a frequency  band  centered  about  15.5  MHz  was  chosen  as 
having  the  best  propagation  from  the  weather  watch  area.  A clear- 
channel  search  was  then  conducted  and  a frequency  of  16.09  MHz  was 
selected  for  use  in  data  collection. 


Appendix  B 

THE  .ASSOCIATION  OF  ATMOSPHERICS  RECEIVED  IN  CALIFORNIA 
WITH  THOSE  RECORDED  UNDER  A NEW  MEXICO  THUNDERSTORM 

In  Section  IV,  atmospherics  received  in  California  were  associated 
with  those  received  close  to  a thunderstorm  in  New  Mexico.  These  two 
sets  of  data  were  recorded  using  different  time  standards.  While  an 
effort  was  made  to  synchronize  these  standards,  it  was  recognized  that 
a relati ve-time  uncertainty  was  unavoidable.  To  remove  this  uncer- 
tainty, identifiable  signatures  of  a strong  atmospheric  received  in 
California  were  matched  to  corresponding  features  of  the  ELF  waveform 
received  in  New  Mexico.  With  this  procedure,  the  time  axis  of  the  New 
Mexico  data  was  shifted  by  11  ms  and  this  shift  was  then  used  for  com- 
paring the  remainder  of  the  five  minute  data  sets.  This  shift  was  made 
in  order  to  compensate  for  the  4.5-ms  difference  in  propagation  paths 
and  for  the  difficulty  of  synchronizing  to  the  weak  and  fading  skywave 
time  signal  received  from  the  National  Bureau  of  Standards. 

The  introduction  of  this  time  shift  introduces  the  possibility  of 
bias  to  the  association  process,  and  it  is  important  to  establish  that 
subsequent  associations  were  due  to  a definite  cause-and-effect  rela- 
tionship and  are  not  a mere  happenstance  brought  about  by  the  arbitrary 
initial  time  shift  and  the  chance  sequence  of  random  atmospherics. 

The  atmospherics  received  in  California  were  presumably  caused  by 
in-beam  lightning  flashes.  Only  a fraction  of  these  atmospherics  were 
likely  to  be  caused  by  the  storm  monitored  in  New  Mexico.  Clearly, 
the  more  atmospherics  received  in  California  and/or  the  looser  the 
rules  for  association  (i.e.,  the  association  time  interval),  the 
greater  was  the  chance  of  a random  association.  Thus  it  was  necessary 
to  determine  the  likely  rate  of  random  associations  based  on  the  rati 
of  atmospherics  received  in  California  and  the  time  interval  used  lor 
associations . 

During  the  5-minute  data  sequence  the  average  rate  of  at  im  > - ; 1 
received  at  WARF  was  less  than  five  per  second.  Here  each  di  t 
atmospheric  signature  that  could  be  used  for  making  an  ass, 

S 7 
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(principally  the  abrupt  return-stroke  quenching)  was  counted  separately. 
While  the  average  association  showed  much  less  than  a 5-ms  discrepancy, 
for  the  purpose  of  this  analysis  a ±0  ms  association  interval  was 
chosen.  Based  on  these  numbers,  the  probability,  p,  of  randomly  asso- 
ciating a New  Mexico  event  wj th  an  unrelated  atmospheric  received  in 


California  is  equal  to  the  rate  of  atmospherics  received,  multiplied 

-2 


by  the  association  interval  (i.e.,  5 x 1 Y 10  0.05).  Assuming 

independent  events,  the  number  of  associations,  k,  out  of  the  total 
number  of  New  Mexico  events,  N,  is  given  in  accordance  with  the  binomial 
distribution : 


r.  (K  associations^  N!  K ,N-K  , . 

p Uut  of  N trie4  = KfTinnoT  P (1  - p)  ‘ (B*1) 


During  the  5-minute  data  set  analyzed,  53  flashes  occurred  in  the 
New  Mexico  thunderstorm,  excluding  the  flash  used  to  adjust  the  time 
scales.  If  none  of  the  atmospherics  received  in  California  were  caused 
by  these  New  Mexico  flashes,  then  Eq . (B-l)  would  describe  the  proba- 
bility of  K random  associations  with  p = 0.05  and  N = 53.  The  proba- 
bility of  having  more  than  seven  random  associations  is  found  to  be 
less  than  0.01.  There  were  a total  of  43  observed  associations.  Since 
there  is  much  less  than  one  chance  in  one  hundred  that  most  of  the  ob- 
served associations  could  have  occurred  by  chance,  the  desired  causc- 
and-effect  relationship  has  been  established. 
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Appendix  C 

THE  WARF  ANTENNA  DIGITAL  BEAMFORMER 


Section  IV  of  the  main  test  describes  observations  of  the  spatial 
distribution  of  individual  atmospherics.  Those  observations  were  ob- 
tained by  digitally  forming  eight  contiguous  beams,  each  of  approxi- 
mately 0.5°  beamwidth.  The  equipment  and  processing  required  to  ac- 
complish this  digital  beamforming  are  described  here. 

Figure  C-l  shows  diagramatically  the  WARF  receiving  antenna  and 
the  equipment  configuration  used  by  the  digital  beamformer.  Elements 
within  each  subarray  are  combined  via  analog  means.  This  enables  the 
array  to  be  steered  in  the  general  direction  of  interest.  The  subar- 
ray beam  pattern  is  also  useful  for  rejecting  unwanted  off-axis  energy. 

43 

The  eight-channel  coherent  receiver  and  simultaneous  sample  and  holds 
were  built  at  SRI  expressly  for  such  digital  beamforming  applications. 
Table  C-l  lists  the  pertinent  parameters  used  for  data  collection  and 
processing . 

The  eight  channels  of  data  were  digitized  and  recorded  for  off- 
line analyses.  During  analysis  the  data  are  read  back  from  tape  and 
the  eight  beams  are  formed  by  applying  a two-dimensional  Fourier  trans- 
form to  the  Data.  Figure  C-2  indicates  the  processing  steps  involved. 
Eight  contiguous  samples  from  each  of  the  eight  receiver  channels  are 
processed  coherently  to  form  a single  azimuth  profile.  A 2048-Hz 
sample  rate  enables  the  calculation  of  independent  azimuth  profiles 
for  every  4 ms  of  data. 

Each  block  of  data  is  first  weighted  to  reduce  sidelobes  in  the 
time/frequency  and  space/azimuth  dimensions.  The  eight  weighted  sam- 
ples from  each  channel  are  then  Fourier  transformed.  This  results  in 
four  complex  frequency  samples  for  each  of  the  eight  channels.  Next, 
for  the  two  center  frequencies,  a Fourier  transform  is  applied  to  the 
data  from  the  eight  channels.  Just  prior  to  this  transform  the  data 
are  padded  with  zeros  to  achieve  azimuthal  interpolation.  After  the 
transforms  the  magnitudes  of  the  data  are  computed  and  the  results 
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ELEMENTS 


Table  C-l 

SUMMARY  OF  DIGITAL  BEAMFORMING  PARAMETERS 


658D32  967 


FIGURE  C-2  BEAMFORMER  DATA  PROCESSING  FOR  EACH  4-ms  DATA  INTERVAL 


from  the  two  frequency  cells  are  incoherently  averaged.  This  procedure 
results  in  an  interpolated  estimate  of  the  azimuthal  distribution  of 
noise  received  during  a 4-ms  period.  Contiguous  4-ms  data  intervals 
are  processed  in  an  identical  manner  to  produce  the  plots  of  amplitude 
versus  azimuth  and  time  shown  in  Figures  15,  16,  and  17  of  the  main 
text . 
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